Approach to a rational rotation number 
in a piecewise isometric system 

J. H. Lowenstein and F. Vivaldif 

Dept. of Physics, New York University, 2 Washington Place, New York, NY 10003, USA 
'\ School of Mathematical Sciences, Queen Mary, University of London, London El 4NS, UK 

Abstract 

We study a parametric family of piecewise rotations of the torus, in the limit in 
which the rotation number approaches the rational value 1/4. There is a region of 
positive measure where the discontinuity set becomes dense in the limit; we prove that 
in this region the area occupied by stable periodic orbits remains positive. The main 
device is the construction of an induced map on a domain with vanishing measure; 
this map is the product of two involutions, and each involution preserves all its atoms. 
Dynamically, the composition of these involutions represents linking together two sec- 
tor maps; this dynamical system features an orderly array of stable periodic orbits 
having a smooth parameter dependence, plus irregular contributions which become 
negligible in the limit. 
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1 Introduction 



This paper is devoted to the study of the action of the matrix 

C={\ ~o^) A = 2cos(2^p) (1) 

on the unit square, in the limit A — > 0. This deceptively simple dynamical system has 
a surprisingly complicated behaviour. In figure [2] we display the phase portrait in the 
North-East corner of the unit square, for A = . The picture suggests the existence of 
an infinite hierarchy of stable islands, immersed in a 'pseudo-chaotic' sea. 

More precisely, we let = [0, 1)^ and we define the map 

{x,y) ^ {\x -y + i{x,y),x) i{x,y) = - [\x - y\, (2) 



Figure 1: Left: Detail of the phase portrait of the map F for A = 2~^, corresponding to an 
irrational value of the rotation number p, close to 1/4. A portion of the large central island 
appears in the SW corner of the picture. Right: Magnified view of the region near the NE corner. 

where [-J denotes the floor function. This map is area-preserving and has time-reversal 
symmetry [24] 

F-i=GoFoG-i G:Q^n {x,y)^{y,x). (3) 

Occasionally, we will regard -F as a map of the torus, rather than of the square. It must 
be understood however, that the function i in (l2|) is defined over Q., not over M^/Z^. 

In the specified parameter range, the map F is non-ergodic, and has zero topological 
entropy [8]. It is linearly conjugate to a piecewise rotation on a rhombus with rotation 
number p, and indeed it's a piecewise isometry with respect to the metric induced by the 
inner product 

Q{U, V) = U,V, + UyVy - ^ {UxVy + UyV^) . (4) 

The quantity Q(C/, U) is the invariant quadratic form of the matrix C in ([T|). For A = 0, 
this metric reduces to the ordinary Euclidean metric. Throughout this paper, the term 
isometry will always refer to the metric Piecewise isometries are dynamical systems 
that generalize to higher dimensions the construct of interval exchange maps. Recently, 
these systems have attracted a great deal of attention, both theoretically [1,3,7,13-15,18, 
20,23,32], and in applications [4,9-11,25,28,33]. 

In the analysis of these systems, a decisive simplifying factor is the presence of self- 
similarity, which allows a satisfactory — occasionally complete — characterization of the 
dynamics [30]. Self-similarity has been invariably found in maps with rational rotation 
number (the parameter p in ([I])), although the occurrence of scaling in this case has never 
been justified by a general theory. For rational rotations, the stable regions in phase 
space — the ellipses of figure [1] — become convex polygons, and the system parameters 
(such as the quantity A in ([T])) are algebraic numbers. Much of recent research has been 
devoted to quadratic parameter values, plus some scattered results for cubic parameters 
[18,21,26,27,30]. Rational rotation numbers with prime denominator were considered 
in [17] from a ring-theoretic angle, in a rather general setting. In all cases in which 
computations have been performed, the complement of the cells, namely the closure of 
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Figure 2: Left: Mismatch between the action of F on 17 and a rotation by 7r/4, for a parameter 
value close to zero. Shown are the first few images of the boundary of il, which are tangent to the 
central island. Right: Dominant features of the discontinuity set in the North-East sector E. 

the so-called discontinuity set (the dark area in figure [1]), has been found to have zero 
Lebesgue measure. However, no general result has been established in this direction. 

The case of irrational rotations — the generic one — has stubbornly resisted attack, 
largely due to an apparent lack of self-similarity. In 1997, P. Ashwin conjectured that for 
irrational rotations, the complement of the set of elliptic islands has positive Lebesgue 
measure, and this measure depends continuously on the parameter [2]. Subsequently, the 
semi-continuity of the measure at irrational rotation numbers was established rigorously 
[16]; however, the important question of positivity of measure remains unresolved, and 
constitutes the motivation for the present work. Ashwin's conjecture is analogous to the 
well-known conjecture of the positivity of measure of chaotic motions in Hamiltonian 
systems with divided phase space. 

To identify elusive instances of scaling behaviour, in this paper we examine irrational 
rotations in the vicinity of a rational one. For values of the parameter A approaching 
zero, the map F approaches a rotation of the unit square about its centre, by the angle 
tt/2 (figure [21). To a first approximation, the mismatch between the actual and the limit 
dynamics causes the images of the boundary of the square (the discontinuity set) slowly 
to build a regular envelope of the central disk, by line segments. These segment intersect 
transversally, forming an orderly array of pseudo-hyperbolic points [31]. It is not difficult to 
see that in the limit the pseudo-hyperbolic points become dense in the four corner sectors, 
and any stable island present in these sectors must be confined within the meshes of the 
grid, which are of vanishing size. On the other hand, at A = the discontinuity set is just 
the boundary of the unit square, and all this complicated structure disappears. 

This paper is devoted to the study of the limit A — > 0. We have two main results. 
Firstly, we construct an induced mapping that describes in this limit the dynamics of F 
outside the main island. This map has interesting properties, summarized in the following 
theorem. 
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Theorem A. For A positive and sufRciently small, there is a mapping L, induced by F 
on a domain A of area 0(A), which serves as a return map for all orbits outside the main 
island, apart from a set of measure O(A^). The mapping L is the composition of two 
involutions; each involution preserves all its atoms, apart from a set of irregular atoms of 
total measure O(A^). 

Theorem A is the summary of theorems [6] and [TTl which describe the structure of the 
involutions comprising L, and theorem 1231 which estimates the measure of the orbits that 
pass through the so-cahed regular atoms of the map L. This resuh leads us to consider 
the symmetric fixed points of L, which form a distinguished family of stable periodic 
orbits of the map F — the regular sequence. Our second result concerns the area of the 
corresponding stable regions. 

Theorem B. As A — > 0"^, the orbits of F that correspond to the fixed points of L form 
an infinite family of stable symmetric cycles, with the property that the total area of the 
associated islands approaches a positive limit. 

The above theorem will appear in the text as theorems [17] and [22l which also give 
detailed information about these periodic orbits. The 'perturbative' approach, whereby 
irregular contributions of vanishing measure are neglected, is an essential element of our 
analysis, since the complexity of the neglected dynamics is considerable. 

This paper is organised as follows. In section [2] we provide the basic definitions and 
constructs. In section [3] we first define the induced map L, and then provide a geometrical 
proof of the main part of theorem A (theorems [6l and [TT]) . This proof is centred around the 
idea of sector maps, continuous transformations of a neighbourhood of an island, which 
are close to the identity. We will show that the dynamics of F is the result of linking 
together two sector maps (see figures [3] and [4]) . 

In section HI we provide an algebraic proof of theorem \TT\ and also derive explicit 
formulae for the atoms of the return maps. These will be essential for the subsequent 
analysis. In section [5] we prove theorem B. We begin by constructing the regular sequence 
of periodic points, which shape the phase space outside the main island (figured]). These 
orbits are parametrized by pairs of atoms of the two involutions comprising the mapping 
L of theorem A. Then we use the formulae developed earlier to establish some asymptotic 
properties of the atoms of the aforementioned involutions. With these formulae, we prove 
that, as A ^ 0, the regular sequence of stable periodic points of L becomes infinite, and 
the limiting area of the ellipses associated with the regular sequence approaches a positive 
limit (theorem [22]) . Our proof is constructive, and we derive an analytic formula for the 
total area. 

In section [6] we complete the proof of theorem A; we prove that the domain of definition 
of the induced map L — indeed the subset of it constituted by the regular atoms of L — is 
a surface of section for all orbits, except for a set of measure 0(A2) (theorem [23]). Finally, 
in section [7] we briefly consider the extension of our results to the case of negative A and to 
the dynamics of irregular atoms. Several formulae and proofs are collected in an appendix. 
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2 Preliminaries 



The map F consists of the action of the matrix ([T]), followed by a translation by l, given 
in equation ([2]). The quantity l takes a finite set of integer values. In particular 

L{x,y)£ {0,1} if 0<A<1. (5) 

We shall now assume that A belongs to this range, approaching zero from above; later on, 
we shall restrict the parameter further, to a smaller right neighbourhood of zero. The case 
of negative parameters will be considered briefly in section [71 

The unit square is partitioned into two atoms ilj, which are convex polygons given 
by the level sets of the function l 

^i = t-\j) j = 0,l. 

Given a cod^ l = {lq, li, . . .), with integer symbols it taken from the alphabet ([5]), we 
consider the set C{t) of the points z £ il. for which l{F^{x, y)) = tt, for t = 0,1, . . .. These 
are the points whose images visit the atoms in the order specified by the code. If C(t) is 
non-empty, then it is called a cell. Generically, a cell is an open ellipse together with a 
subset of its boundary [16, proposition 2]. In coordinates relative to their centre, these 
ellipses are similar to the ellipse Q{x,y) = 1, cf. equation (jl]). A cell can also be a convex 
polygon (when the rotation number p is rational), or a point. 

The set of all images and pre-images of the boundary of the atoms constitutes the 
discontinuity set T 

oo 

r= y F\dVL) dn = [jdni 

t= — OD i 

which consists of a countable set of segments. For the purpose of generating T, the set 
d^l may be replaced with any set that covers di} under iteration. One verifies that the 
discontinuity set of the map ([2]) is generated by the oriented segment 7, with first end- 
point (0,0) (included) and second end-point (0,1) (excluded). The map ([2]), viewed as a 
map of the torus, is discontinuous only on 7, from the left. Indeed 

F(l -x,y)- F{x, y) = (A + i(l - x,y) - i{x, y), 1 - 2x) = (A, -2x) (mod I?). 

Letting x — > 0"*", we see that F has a sliding singularity at 7 [22], which causes a shift 
by A in the horizontal direction. By time-reversal symmetry, the map F^^ has a siding 
singularity on (^(7). 

The generator 7 is oriented in such a way that the map F be continuous on the 
segment's right side. The images and pre-images of 7 inherit an orientation with the same 
property. In view of this, we'll say that a segment in the discontinuity set is glued to any 
domain tangent to it on the right. Furthermore, if the first and second end-points of a 

^We use the symbol t for both the coding function and the code. 
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segment belong to two sets A and B, respectively, we say that such a segment connects 
A to B. 

Necessarily, any periodic orbit of F will have a point whose cell is tangent to 7, from 
either the left or the right. For t G Z, the set F^{'y) is a collection of segments. Given a cell 
£, the regular component of -F*(7) with respect to £ is the set of segments in -F*(7) that 
are tangent to £. A non-empty regular component consists of a single segment; indeed 
all components of the image of 7 must have the same orientation, so that two tangent 
segments would have to lie at opposite points of £. However, for A 7^ 0, there is no cell 
touching 7 on opposite sides, and the image of a tangent point is a single point. The same 
property holds also for images of the segment G{'y), and indeed of any other segment. We 
will also speak of the regular component with respect to a cycle, because, given a cycle 
and an integer t, there is at most one ellipse in the cycle such that a component of -F*(7) 
is tangent to it. Regular components will be important in section [3l in the construction 
of the induced map. 

At A = the dynamics is trivial. The atom 0,q (a triangle) collapses to the segment 
G(7), and the discontinuity set is just F = 7U G^j). All points in F have period 2, except 
for a fixed point at the origin. The rest of the space — the interior of O — consists of a 
single square cell; its points have period four, except for a fixed point at the centre. 

We are interested in periodic orbits for parameter values near zero. A periodic orbit 
(xo, xi, . . . , xt-i) corresponds to a periodic code i = {lq,li, . . . , it-i)- As the parameter 
varies, the points of the orbit move in phase space; in [31, theorem 2] it was shown that 
the coordinates Xj are rational functions of A. Specifically, one has 

-^(^^^) = ^^xm^ ^■=0,.. (6) 

where Xt and Ait are polynomials in Z[A], and a is the left shift map. The denominator 
Mt is monic of degree [{t + 2)/2j and it depends only on the period; the numerator Xt 
has degree [{t + l)/2j , and it depends on the orbit via the code i. This algebraic structure 
is common for periodic points of piecewise affine systems, see, e.g., [5,6]. In section [5] we 
will derive ad hoc expressions for these functions, tailored for the A ^ regime. 

The symbolic dynamics is far from being complete, so not all periodic codes l corre- 
spond to an actual periodic orbit. For this to be the case, the value of all functions Q 
must belong to the half-open interval [0, 1). The function that is closest to the boundary 
of the unit interval determines the common size of all cells of the orbit corresponding to 
that code. Accordingly, we define the radius r of a periodic orbit as 

r(A) = min {x,(A), 1 - x,(A)}. (7) 

o<i<t-i j\ /J 

The function r is piecewise rational; besides the singularities inherited from the functions 
Xj, the radius will typically have points with discontinuous first derivative, due to a change 
of the index j for which the minimum in ([7]) is attained. An orbit exists if it has positive 
radius, and the radius becomes zero at the bifurcation parameter values. An orbit of zero 
radius exists if and only if none of its rational functions assume the value 1, which is the 
excluded point in the unit interval. On the torus, this distinction becomes irrelevant. 
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Formulae of the type ([6]) also describe the parameter dependence of pseudo-hyperbolic 
points [31, theorem 6]. These are the points that recur to the boundary of the atoms. 
More precisely, a point in Q is pseudo-hyperbolic if it maps to 7 in the forward time 
direction, and to G{-f), in the backward time direction. For an irrational rotation number, 
these points correspond to transversal intersections of two segments of the discontinuity 
set, which act like 'pseudo-separatrices'. It can be shown that, generically, an orbit on 
the discontinuity set can contain only finitely many pseudo-hyperbolic points, the first of 
which is in G{'j), and the last in 7. A pseudo-hyperbolic point can be periodic only for 
parameter values which are algebraic of degree greater than one [31, theorem 9]. Thus, 
in general, a pseudo-hyperbolic sequence belongs to an infinite orbit; such a sequence 
is determined by a finite symbolic dynamics, which describes the itinerary between the 
endpoints on G{'y) and 7. 

Pseudo-hyperbolic points will occur as vertices of the atoms of the return map L, to 
be defined in the next section. As with periodic orbits, we shall derive explicit formulae 
for these points, appropriate for the A — > limit. 



3 Return map 

In this section we prove the main part of theorem A. We will construct an induced mapping 
L that describes the dynamics of F in the limit A ^ 0^. This map is the composition 
of two involutions on a A-dependent domain whose area vanishes in the limit (theorems [6] 
andE]). 

From the time-reversal symmetry ([3]), we find that F can be written as the composition 
of two orientation-reversing involutions 

F = HoG with H = FoG, G^ = = Id (8) 

with H given by {x,y) 1— > ({Ay — x}, y), where {•} denotes the fractional part. Let FixG 
and FixH be the sets of fixed points of the involutions G and H, respectively. The set 
FixG is the segment x = y, independent of A; for < A < 1, the set FixH consists of two 
segments, with end-points (0,0), (A/2,1), and (1/2,0), ((A + l)/2, 1), respectively. 

We begin by examining the rational functions of the main periodic orbits of F (cf. equa- 
tion dll)) 

t code denominator numerators radius 

1 (T) 2- A (1) l-xo(A) (9) 

2 (01) (2, A) XI (A) 

These cycles are symmetric (G-invariant), and the 2-cycle belongs to FixH. At A = 
the fixed point has radius 1/2, while the 2-cycle has zero radius, but it still exists on the 
discontinuity set: {(1/2,0), (0, 1/2)}. It turns out that, in the A ^ 0^ limit, the dynamics 
is dominated by these two periodic orbits, which determine the coarse features of the 
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phase space. The fine structure also depends on other periodic orbits, some of which will 
be analysed in section [71 

The ellipse E bounding the large cell of the 1-cycle is tangent to the segment 7 at the 
point (1,t), where r = (1 + A)/2. By G-symmetry, E is also tangent to 70 = G(7) at 
the point Tq = (r, 1). The open region delimited by the segments 7, 70, and the elliptical 
arc connecting Tq to G(7b) clockwise will be called the 1-sector, denoted by S. By 
construction, the set S is invariant under G. The 1-sector — indeed a vanishingly small 
subset of it — will serve as a surface of section for orbits outside the main island. Not all 
such orbits intersect E, for example, the 2-cycle displayed in Q, and the cycles mentioned 
in section [71 However, the missing orbits will be shown to have vanishing measure. More 
important, the distinction between the orbits that enter S and those that don't will be 
essential to our analysis. 

For A near zero, the map F'^ is close to the identity, and hence the set S is 'close to' 
being invariant under F^. There are however orbits that enter and leave S; to characterize 
them, we define the sets 

A = S\F^(S) H = F^(S)\S (10) 

where, by time-reversal symmetry, 

G(A) = S\F-^(S) H = F^oG(A). 

Thus G'(A) is the set of points that leave S under F'*. Below, we will define two induced 
maps over the set A. 

To determine the set A, we construct the first few images of the segment 7, which are 
listed in table [H below as pairs of end-points. We adopt the following conventions: 1) an 
integer n in boldface denotes the n-th iterate of the generator 7 (e.g., 3 = ^'^(7)); ii) 
segments crossing 7 have negative x-coordinate in the end-point on the left of 7; iii) if an 
end-point of n belongs to more than one set m, we list only the value of m corresponding 
to the lowest order iterate. 



segment 


first 

end-point 


second 
end-point 





(0,0) 


(0,1) 


1 


(1,0) G 


(0,0) G 


2 


(A,l) G 1 


(0,0) G 


3 


(A2,A)g2 


(1,0) GO 


4 


(-A + a3,a2) e 3 


(A,1)G1 




(A-2A2 + A^1- A + A3) g4 


(5-1,1) G 1 




(<?',o)gi 


(A2,A)g2 




1 + 2A - A2 - A^ , 1 


-A2 



^ 2-A2 f 2-A2- 
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Figure 3: The sector E, and its image under F'^, corresponding to a small clockwise rotation. 
The darker region is the entry sub-domain A; the exit sub-domain G(A) is the reflection of A with 
respect to the symmetry axis Fix G (the dashed line) . All points in E will eventually map to the 
region E = F'^ o G(A), located on the opposite side of the discontinuity hue 0. 



Table 1: The first five iterates of the generator 7 of the discontinuity set. 

The segments 0, . . . , 4 are regular components with respect to both the 1-cycle and 
the 2-cycle. Since < A < 1, the segment 4 crosses the discontinuity hne 7, and so its 
image consists of two segments 5^ and 5^. These are the images of the portions of 4 lying 
to the left and to the right of 7, respectively, and they are regular components of 5 with 
respect to the 1-cycle and the 2-cycle, respectively. 

We define approximate orthogonality relations. If the angle between two segments m 
and n is equal to 0, tt -|- 0(A), we say that m and n are quasi-parallel, and write m || n. 
Likewise, two segments are quasi-perpendicular (m _L n) if the angle between them is 
ib7r/2 -|- 0(A). Plainly, m and n are quasi-parallel (perpendicular) if m — n is even (odd). 

There are also incidence relations. We use the notation m H n to indicate that one 
endpoint of m lies in the interior of n (and similarly for m h n) . 

The images of 7 have the following behaviour under symmetry 

0(m) = GoF™(7) = F-"oG(7) =F-"oF(7) = -m+l (12) 
H{m) = FoF™(7) = Foi7oGoF™~i(7) = GoF"'-i(7) = -m + 2. (13) 
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We specify polygons/sectors by listing their boundary components between angle brackets, 
arranged clockwise, e.g., 

Oo = (0,l,-1) S = (f,l,0). 

Since G and H are orientation-reversing, mapping a polygon entails mapping its sides, 
and reversing their order, e.g., 

G(l^o) = (G(0),G(-1),G(1)) = (1,2,0). 
The segment 5"*^ connects 4 to 1, and 5^ H S is the segment with end-points 

(1,/), (5,1) 9" = Y^Y^ ^^^^ 

where g was given above. As A approaches zero, 5 — > 1/2 and g" 1. For our purpose it 
is sufficient to require that g" < g; so we replace the bounds in ([5]) by the more restrictive 
conditions 

0<A<A+<1 where A+ = 2 cos(27r/9). (15) 
(The number A+ is algebraic, being a root of the polynomial x"^ — 3x + 1.) 

The segment 5^ is tangent to the ellipse £ at the point Ti; hence 5^ decomposes S into 
three regions (figure [3]); a large triangular sector F^(S) n S, of area 0(1), a right triangle 
of area 0(A), and a small triangular sector e of area 0(A3). The set A is the union of the 
last two regions. We find 

A=(0,5\l)Ue H = (0,5\4). (16) 

Our analysis will be perturbative, and in what follows we shall neglect domains of area 
O(A^). In particular, we shall omit e from all considerations, and represent A as a triangle. 

We intend to study the dynamics outside the main island via the first return map L 
induced by F on A. This map is constructed as the composition of two transit maps 

L = Z°"* o Z''^ Z'"^ : A ^ G(A) Z"*^* : G(A) ^ A (17) 

where Z™ and Z"*^* are the first-return maps induced by F on the respective sets. The 
map Z"^ is built up from iterations of F'^ inside S, the map Z°"* is built up (after the 
initial step out of S) by the iteration of F outside S. 

The definition of the maps Z™'"*^* requires some care. As noted above, some orbits 
outside the main island do not enter S at all, and it is conceivable that some orbits exiting 
S will never return there. From Poincare recurrence theorem, we know that these points 
have zero measure; so we'll ignore their contribution and accept that the map Z""* may 
be undefined on a zero measure subset of G{A). In particular, the transit time of the map 
L°^^ may be unbounded. The same is true for Z™, as a result of having dropped the small 
set e from the domain A. 

The following result has important consequences for the dynamics of the sector S. 
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Lemma 1 The restriction to the sector S of the maps F^^ is continuous. Continuity 
extends to the boundary of apart from a sliding singularity of on 7, and of F~^ on 
G{'y). The symbolic dynamics is (1, 1, 1, 1) for points in S \ G(A) and (1, 1, 1, 0) for points 
in G(A). 

Proof. Since F is continuous away from 7, and 7 = is a sliding singularity (the 
segment does not detach itself from S under iteration) we have 

= (£:,t + l,t), t = 0,...,3 F^{^) = {£,5\4). (18) 

From table [H we find that, in the parameter range (jlSp . of the first four images of 0, 
only 4 intersect 7. However, the point at which 4 is tangent to £ lies to the left of 7. 
Hence none of the first three images of S intersects 7, and therefore F^ is continuous on 
S. The continuity of the inverse follows from time-reversal symmetry. Prom (jlSp . we find 
that -F(S) n S = 0, and hence all points in S \ G{A) return to S in four iterations of 
F and no fewer, with return symbolic dynamics (1,1,1,1). Furthermore, from (jl3p and 
(jl6l) . we have G(A) = (—4, 1, 0), and since G(A) C S and F^ is continuous there, we have 
F3(G(A)) = (-1,4, 3) C r^o- So the symbolic dynamics of a point in G{A) is (1,1,1,0). □ 

The transit maps are piecewise isometrics, and our next step is to partition 

their respective domains A and G(A) into atoms. Since the two domains are mirror image 
of each other under the involution G (see equation (jlOp ). we replace ()17p by the more 
convenient decomposition 

L = L°"* o V""; L'°'°"* : A ^ A, = G o L'°, = L°^^' o G (19) 

where A is the common domain of the two maps. We will often represent the right triangle 
A using scaled coordinates, with the shorter leg being four times the length of the longest 
one — see figure [71 

The distinction between regular and irregular atoms will be important in what follows. 
All regular atoms share a simple geometrical structure (apart from two exceptions), de- 
termined by regular components with respect to a periodic island. Asymptotically, these 
atoms occupy the entire domain A™ of definition of the maps. Irregular atoms are more 
varied, but their detailed description is not required in asymptotic calculations, because 
in the limit, their total area becomes negligible. 

3.1 Preliminary lemmas 

The proof of theorem A will be quite laborious. We begin with some preliminary lemmas 
on involutivity and symmetry. We call an atom maximal if it is not properly contained 
in another atom. 

Lemma 2 Let L be a piecewise isometric involution. Then the image of a maximal atom 
is an atom. In particular, a maximal atom containing a fixed point is invariant. 
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Proof. Since L is an involution, its atoms are the same as those of L^^. If $7 is an 
atom of L, we claim that L is continuous on L{Q). Indeed if it were not, we could find 
arbitrarily close points z,w £ L{Q) mapping far apart in thereby contradicting the fact 
that L is an isometry on 0,. Thus L(0,) C O', where J7' is another atom. If this inclusion 
were proper, then L{Q') would be a connected set properly containing 0, because L is 
continuous on f2' and volume-preserving. But since L is continuous on L{i}'), then 
would not be maximal. So, if is maximal, then = fi', as desired. Furthermore, if 

contains a fixed point of L, then Q intersects L(0,), and is therefore invariant. D 



Lemma 3 The maps L™'"*^* are involutions. 



Proof. Let z £ A. Since is an induced map, we have Z™(z) = F^(z), for some 
t = t{z). From ([3]), a straightforward induction gives = G o o G, for all t £ 1^, and 
hence 

L"^or°(z) = GoF*oGoF*(z) =F-*oF*(z) = z. 

Since z was arbitrary, is an involution. The involutive nature of is proved with 
an analogous argument. D 



Next we look at the fixed sets of these involutions. 



Lemma 4 Let A™'"*^* be an atom o/L™'°"*, and let t = t{k) be its transit time (number 
of iterations of F). Then 

Yh^Lf = F-'{¥\^G)r\Kl if t = 2s 

FixLjf = F-' oG{Yi^H)r\Kl' if t = 2s + l 

FixL^"*= F"(FixG) n A^"* if t = 2s 

Fix Ll''^ = F' (Fix H)nAl''^ if t = 2s + 1 

where ^5 the restriction of to A™'""'. 



Proof. From p9|) . the fixed points equations become 

L^{z) = z^Go F\z) = z L^'^'iz) = z^F^o G{z) = z (20) 

and we are interested in solutions (if any) within A™'°"*. First, we consider the inner map 
L™, and z G A'^'^. If i = 2s, equation ([20|) becomes 

GoF^'{z) = z <^ F-^'oG{z)=z 

4^ F-' o G{z) = F%z) 
^ GoF^z) = F'{z). 
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This shows that F*(z) G FixG, namely that z G F~''(FixG). If t = 2s + 1 is odd, we 
obtain, using ([8]) 

^ FoF'{z) = GoF'{z) 

4^ HoGoF\z)=GoF'{z). 

We find that G o G Fixff, or z G F"* o G(Fixi7). 

Similarly, for z G A^"*, and t even, the outer map gives 

F^'oG{z) = z ^ F-' o G{z) = F-%z) 
<^ GoF-%z) = F-'{z), 

so F-^{z) G FixG, that is, z G F*(FixG). Finally, if t is odd 

F^'+^ o = z ^ F'+^ o G(z) = F"'{z) 

^ GoF-^oF-'{z)=F-'{z) 
^ HoF-%z) = F~'{z) 

which gives F^''{z) G Fixff, that is, z G F'*(FixF). □ 

We will show that for a significant class of atoms, the transit times of and 
are even and odd, respectively, and the fixed lines of the involutions intersect these atoms, 
which are themselves invariant. This will be done in the next two sections, where we look 
in detail at the transit maps. 



3.2 Structure of U 



The map results from restricting to S the action of the matrix G^ (see equation ([TJ) 
with respect to the fixed point of F. We consider again the point Tq at which the ellipse 
£ is tangent (on the torus) to the segment 1, see figure [3l Let Ti, T2, . . . be the images of 
To under F^, and let 



N 



1 



1 



4(1 - 4p) 4 



— cos-i(A/2), 



(21) 



where p is the rotation number, see equation ([T]). Then the first N points of this sequence 
lie on the ellipse f , within the boundary of S. 



For n = 0, ... ,A'^ — 1, we let 7„ = (4n + 1)^ be the regular component of F^"(l) 
with respect to £. Then 7„ is tangent to £ at the point T„. For any k, 7„ and ^n+k are 
quasi-parallel, and their orientations agree. The segments 7„ will be used to construct the 
atoms of L™. We begin with a lemma, whose geometric content is illustrated in figure [5j 
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Figure 4: The linked sector maps of the 1-cycle (the large circle) and the 2-cycle (the small circle), 
which lie on opposite sides of the discontinuity line 0. Under the action of the map F**, part of 
the sectors E and E' (the darker regions) escape to the triangles S and S', respectively, located on 
the opposite side of 0. The points in S that do not end up in E' (under F^) comprise the triangle 
Ti; those that do comprise the quadrilateral T2. 

Lemma 5 For n = \, . . . ,N — \, the regular component 7„ connects 4 to jn-i, o-nd it 
intersects transversally the segments and —4^ on the boundary of G{K), apart from a 
possible non-generic degeneracy o/7jv_i n —4^ into a segment. 



Proof. From (jl6p we verify that and —4^ do indeed belong to the boundary of G(A). 
Let n = 1. The segment 71 = 5^ connects 4 to 70 = 1 — see table [TJ For the parameter 
A in the range ([IS]), the point Ti lies outside G(A), and therefore 71 crosses and —4^ 
transversally. We proceed by induction on n, and assume that for some n in the range 
1 < n < A'^ — 1, the segment 7„ connects 4 to 7n-i, and it is tangent to f at a point 
Tn lying outside G(A). Then 7„ intersects 7 transversally, so we consider the segment 
7~ = 7n n S. By construction, the segment 7" connects 7 to 7n-i- From lemma [H the 
set 7n+i = F^iln) is a segment tangent to £^ at a point T„+i, which connects 4 to 7^ (the 
statement concerning the first end-point is justified by the sliding nature of the singularity 
on 7). Since is close to the identity, the orientations of 7„ and 7^+1 agree. Because 
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G(A) 




Figure 5: The geometry of lemma [5l the regular components 7„ intersect transversally the bound- 
-4^ and of the exit domain G(A), creating pseudo-hyperbolic points. 



aries 



71 < A^ — 1, the point T^+i lies outside G'(A), which completes the induction. An exception 
may occur for n = — 1, in the non-generic case in which Ttv-i hes on the boundary 
segment —4? of G(A). In this case 7Ar_i is tangent to that boundary. □ 



Prom the above lemma, we see that the segments 7„ decompose G(A) into — 1 
quadrilaterals, plus a small residual region near the tip. We are then led to consider the 
half-open quadrilaterals A^ C A, given by 

Ajf = (70, G(7„-i),7i,G(7n)) 

= (l,(-4n + 4)\5\(-4n)i) n = 1, . . . , A^ - 1. (22) 

The sides are oriented in such a way that 

70, G(7„-i) n Aj,'^ = 71, G(7n) n A^'^ ^ 0. (23) 

We also represent polygons as lists of vertices 

AL" = [Qo(n),go(n-l),Qi(n-l),Qi(n)] Qo,i(n) = 70,1 n G(7n). (24) 

Below, we will show that the sets A™ are atoms of L™. Since these domains are bounded 
by regular components, we shall call them regular atoms; all other atoms will be termed 
irregular. 

From ()24p . we see that the vertices of A™ are pseudo-hyperbolic points of the map F. 
Indeed, forra = 1, . . . , A^— 1, the point Qo{n) is the intersection of 70 = ^(7) and a segment 
in the discontinuity set, and hence it is the left end-point of a pseudo-hyperbolic sequence 
(see remarks at the end of section [2]). The point Qi{n — 1) = F'^{QQ{n) belongs to the 
same orbit, and since it is the transversal intersection of two segments (by virtue of lemma 
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[5]), it also is pseudo- hyperbolic. Under the action of F^, these points propagate across 
the sector S, forming a regular pattern of intersecting pseudo-separatrices — see figure El 
The point G{Qo{n)), which belongs to 7, is the right endpoint of the pseudo- hyperbolic 
sequence that starts at Qo{n). 

Theorem 6 The regular atoms of the map are the half-open quadrilaterals \22^) . with 
the boundary specified in [23\). For each n, the atom A™ is invariant under L™, with return 
symbolic dynamics (the exponent denotes repetition of symbols). The fixed set 

of is the segment connecting the vertex Qi{n) to the opposite vertex Qo{n — 1), where 
<5o,i is defined in {2^^ . There are at most five other atoms, of total area 0{}?). 



Proof. Consider the quadrilateral 

G(Ajf) = (7n-i, G(7o), 7n, ^(71)). 

From the proof of lemma O we know that 7„ and G(7i) are glued to G(A™), but the 
other two sides are not. Applying _F~^("-~i) to each segment in G(A™) we obtain the 
quadrilateral (70, G(7„_i), 71, G(7„)), and since the iterates of these segments remain 
within S, we deduce from ()22p that A™ contains an atom of L™, glued to 71 and G(7„), 
with transit time 4(n — 1). 

It remains to show that A™ is actually an atom. The case n = 1 is trivial, since 
A™ = A n G(A), with zero transit time. For n > 1, we must verify that, in the recursive 
construction of A™, no forward image of the segment 7^_i = 7n-i H H (the piece of 7n-i 
that gets cut off by the discontinuity) reaches G(A™) before 7„, namely that 

4 

Ui^*(7+-i)nG(Ajr) = n>l. (25) 
t=i 

Let z E '^n-v From lemma [5] and tabled! we have z E G{Q.q). We have two cases. If 
F{z) E $^0) then F'^{z) E G{Q.q), and F^{z) is either again in VLq, or in ( —1, 2,3); in either 
caseF'^(z) is outside E, hence outside G(Ajf). IfF(z) E ili,thenF(2) E ( -1, 2, 3 ), hence 
F'^{z) E (0,3,4). The points of the latter triangle map either to G(i7o) or to A. From 
neither set it is possible to reach G(A™) in one iteration. We have established equation 
([25]) . and A™ is an atom of L™. 

Generically, the point T/v lies between G{Tq) and G(Ti), and the residual region con- 
sists of one quadrilateral, two triangles and two triangular sectors. At isolated parameter 
values (determined by the condition that 1/4(1 — 4p) be an integer), we have Tn = G{Tq), 
and the residual region reduces to one triangular sector, while the regular atom A^v de- 
generates to a triangle. By construction, the total area of the residual region is O(A^) in 
both cases. 

From ([23]) . we see that 71 and G(7„) are part of the included boundary of A™, and since 
7„ = F'^^"'~^\ji) , the transit time of this atom is equal to 4(n — 1), which is even. With 
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the notation of lemmaHl we have s = 2{n — 1), and hence FixL™ = F"^^^ "^(FixG) D A™. 
Now, by symmetry, 

7nnG(7„) G FixGnS n = 0,...,[(7V-l)/2j. (26) 

We also have the symmetric points 

f4"-2(7o) n G o f4"-2(7o) e FixG n = l,...,[{N- l)/2j (27) 

which he in the South- West sector. (With a shght abuse of notation, we have denoted by 
i^^"~^(7o) the regular component of that set, with respect to £.) Applying the map F~^ 
to the points ([27|) , we obtain a second sequence of points in S 

7n-inG(7„) EF-2(FixG)nS n = 1, . . . , [(A^ - l)/2j . 

Then, by applying repeatedly the identity 

F-4(7,nG(7fc))=7i-inG(7fc+i) 

to pairs of adjacent points in the above two sequences, we translate images of FixG 
within A, and we see that FixL™ is the segment connecting 70 H G(7„_i) to 71 n G(7„). 
Comparison with (j22[) shows that this is the segment connecting the North-East and 
South- West corners of the nth atom (with the obvious modification if the last regular 
atom degenerates into a triangle). 

Now, is an involution (lemma [3j), and its regular atom A™ is maximal, being 
bound by images of the discontinuity line. Then every regular atom is invariant, because 
it contains an invariant segment (lemma [2]). D 



3.3 The sector map of the 2-cycle 

The dynamics of the map L°^^ is dominated by the 2-cycle of F given in table ([9]) , which 
plays a role analogous to the 1-cycle for the map L^°. The 2-cycle generates the regular 
atoms of the map. There are also irregular atoms, but they have negligible measure. 

We noted that the 2-cycle is symmetrical, and it belongs to FixH (see beginning of 
section [3]). We consider the element of the cycle that lies within the atom l^i of F, denoting 
by £' the ellipse bounding its cell (figure H]). From the data ([9]), we find that the segment 
7 = is tangent to £' at the point (0, 1/2). By if-symmetry, the segment -^(7) = ^^(0) = 
2 —see ([13])— is also tangent to £' at the point = H{0, 1/2) = (A/2, 1/2). The open 
region delimited by the segments and 2, and the elliptical arc connecting Tq to H(Tq) 
clockwise will be called the 2-sector T,' 

S' = (<?', 2,0). (28) 

This region is i7-invariant, and both and 2 are are glued to it. 
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The analysis of the dynamics of the 2-sector now proceed as for the 1-sector; we'll use 
the same notation, with primed symbols. First we define entry and exit domains — the 
analogue of ([10]) 



A' = 



F\^')\^' = F^oH{A'). 



Let 



M = 



1 



1 



vr 



as A 



(29) 



2(1 - 4/9) 2 



L2A 







where p is the rotation number; this is the analogue of equation (j2ip . For m = 0, . . . , M—1, 
we let 7^ = (4m + 2)^ be the regular component of F'*™(2) with respect to £' , denoting 
by the corresponding point of tangency. The segments are pairwise quasi-parallel, 
and their orientations agree. The following result is the twin of lemma [H 

Lemma 7 The restriction to the closure of the sector S' of the maps F^^ is continuous. 
The symbolic dynamics is (1,0,1,0) for points in S' \ H{K') and (1,0,1,1) for points in 



Proof. From ^ we find that F(S') = (£:',3, 1). Considering that F(r^) = ((1 + 
A2)/2, A/2) e VLq, from tablefflwe deduce that F(I;') C VLq. Then we have F2(S') C G{VLo), 
hence F^(S') n 7 = 0, and so F^ is continuous on S'. By time-reversal symmetry, the 
inverse F~^ is also continuous on S'. Since the boundaries of S' are glued to it, and 
remain glued under iteration, continuity extends to the boundary. 

Every point in Ti'\H{A') follows the 2-cycle, and so the return dynamics is (1,0,1,0). We 
determine the itinerary olH{K'). We have A' = (0, 6^, 2), and hence H{A') = (2, 0, -4^) e 
VLi. This gives F o H{A') = (3,1,-3) C Oq, whence o iJ(A') C Oi, and finally 
F^ o H{A') = (5^,3, —1) C 0,1. All inclusion relations may be verified with the aid of 
table m together with ((El) and ([H]). □ 

The continuity of F"^ on S' implies the existence of an orderly array of regular compo- 
nents. 

Lemma 8 For m = 1, . . . , M — 1, the regular component 7^ intersects 7 transversally, 
apart from a possible non-generic degeneracy of j'pj_^nH('y[) into a segment. Form > 1, 
7(„ connects 7m_i to 4. 

Proof. The proof of this statement is analogous to that of lemma O with the added 
simplification that continuity now extends to the boundary. The same applies to the 
argument concerning the degeneracy for 'j'j^.j^i. We omit the details for the sake of brevity. 



HiA'). 



□ 



The following result is the analogue of theorem [6] for the 2-sector S'. 
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Lemma 9 For m = 1, . . . , M — 1, the regular atoms of the transit map A' H{A') 
induced by F are the half- open quadrilaterals 

A'^ = (7^,, i^(7™-i), I'l, H[iJ) = (2,(-4m + 4)2,62,(-4m)2). (30) 

The boundary sides 7q and H{'j'j^) belong to A'^, the other don't. The return symbolic 
dynamics of the m-th atom is ((1,0)^™'"'^). There are at most five other atoms, of total 
area O(A^). 



Proof. The proof is analogous to that of theorem [6l with H in place of G, and primed 
variables. Identity (f25|) now reads 

4 

U F*(7„_i n Xi') n A'^ = m > 1. 

t=l 

To establish this identity, we note that 7^_i intersects 7. The right component 7^„x ^ ^' 
will iterate to produce 7^, and we must ensure that the left component 7^_i n H' does 
not interfere with this process. We find 

H'i1^'H'nl^iCF-i(A) = (-l,4,0) 

and the map is continuous over the closure of that domain, from lemma [H apart from 
a sliding singularity on —1 and 0. Thus 

f1(H;)c(0,5\1) f2(H'i)c(1,6\4) 
f3(H'i)c(2,7\3) f4(H'i)c(3,8\4) 

all of which lie outside H', from tabled! Similarly, we find 

r!^ = r.'nno C (-1,0,3) 

and one verifies that F^ is continuous on the closure of that set. Hence 

fHHo) C (0, 1,4) F\Eo) C (1, 2, 52) F\~o) C (2, 3, 6^). 

None of these sets intersects H'; the last one is adjacent to H' an intersects 7. Thus -F^(So) 
is contained in the union of two triangles, which clearly lie outside H'. CI 

As for the sector S, the vertices of the atoms A^ of the transit map over E' are pseudo- 
hyperbolic points of the map F. Indeed such vertices result from transversal intersections 
of regular components — see equation (j30p and figure [H Two vertices of the atom A'^ 
belong to = 7, and hence they are right end-points of pseudo-hyperbolic sequences. 

3.4 Structure of L°"* 

To construct the atoms of the map L°"*, we need to connect dynamically the exit domain 
of the sector S to the entry domain of the sector E', and vice- versa. These connections will 
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not match perfectly, but the mismatch will involve sets of negligible measure. In addition, 
there will be one atom of by-passing T,' altogether (see figure H]) . 

All points leaving S via A"*^* end up in H, and since we want to identify the points that 
will enter T,', we examine inverse images of the entry domain A'. The second pre- image 
of A' intersects H, and so in order to connect the 1-sector to the 2-sector, we consider 
the intersection and symmetric difference of the sets H and F~^(A'). Likewise, connecting 
back the 2-sector to the 1-sector will involve the sets S' and F^^[A). 

Our first lemma quantifies these intersections and symmetric differences. Its proof will 
also provide information about the atoms of L°"*. 

Lemma 10 The points o/A°"* return to E either via the domain Ti = (—2,5-^,4), or via 
A', apart from a set of points of measure O(A^). 

Proof. First, we connect H to A'. By the continuity of F^"^ on S' (lemma [7]), we have 
F-2(A') = (-2,0,4). Defining 

Ti = (-2,5\4) T2 = (-2,4,0,5i) (31) 

one verifies that 

E A F-2(A') = Ti U (-2, 5\0) En F'^iA') = Tg. (32) 

From the orthogonality relations —2 || 0_L5^ || 1, and the fact that the length of the 
side of (-2,5^0) is 0(A), we conclude that |(-2,5^0)| = 0{X^), where | • | denotes 
the two-dimensional area determined by the metric Q. Thus 

\E A F-2(A')| = |Ti| + O(A^). (33) 

We must now consider the dynamics of Ti and of T2, from E to F^^(A). We begin 
with the former. From (jl6p and (j32p . we find, using continuity 

f2(Ti) = (0, 7+, 62) F-^A) = (-1, 0, 4) 

where 7"*" = ^^(5"'^ H E). Hence ^^(Ti) n F^^{A) 7^ 0, and so some points in Ti reach A 
in 3 iterations. However, not all of them do. To see this, we note that 5^ intersects FixH 
transversally, and this intersection belongs to the boundary of Ti. Then, by symmetry, 

Ti Ai7(Ti) = T; Ui7(T;) where = (-2, 5\ -3). (34) 

The i?-symmetric set ()34p will be called a turnstile (figure [6|); by construction, \Ti\ = 
\H(T[)\. Then the set ^^(T'^ U H(T[)) is also a turnstile, which lies on the boundary of 
F-i(A). Finally, 

T'lCTi F\T[) n F-\A) = 
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Figure 6: Left: Symmetrisation of the set Ti (the thick-sided triangle), with respect to the 
symmetry Une FbcH (the dashed Hne). The two darker triangles represent the turnstile. Right: 
Symmetrisation of the quadrilateral T2, with the appearance of two turnstiles; the arrow denotes 
the location of the smaller turnstile. The dashed line is Fix_F^^ o H. The symmetrised polygons 
(light grey) are isometric images of the first and second regular atom of , respectively. 

which identifies T'^ as the set of points of A that map into Ti, but do not return to A 
with it. From the relations —3 || 5^ and the fact that the length of Ti n 5^ is 0(A), we 
have \T[\=0{X^). Thus 

|F3(Ti)nA)| = |Ti| +0(A3). (35) 

There is an analogous phenomenon in the dynamics of T2, defined in equation (j3ip . 
Since H{0) = 2, we have F~'^ o H{0) = —2, and F^'^ o if is an involution. The boundary 
of T2 intersects Fix {F~^oH) in two points, generating two turnstiles (figureEl right). We 
are interested in the intersection 5^ H Fix(F~^ o H), which corresponds to the turnstile 

r'2 U F-^ o H{r'2) where T'2 = (0, 5\ -7). (36) 

Then the set F^(T'2 U F^^ o H{T'2)) is also a turnstile, which lies on the boundary of 
F-i(A). Finally, 

T'2 C T2 F^{T'2) n F-\A) = 

which identifies T2 as the set of points of A that map into T2, hence to A', but do not 
return to A with it. From the relation 5^ || —7, and the fact that the length of T2 n 5"^ 
is O(A^), we have [Tg] = O(A^). This estimate, together with equations ([32]) . and the 
estimates (j33|) and (p5|) . give the desired result. D 

We now identify the regular atoms of Let Ti^2 be as in (fSTI) . and let A^ be as in 
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lemma [H We define the sets 

A°"t = GoF-^(Ti ni7(Ti)) = (0,8,7,1) n A 

AO'it = GoF-4(T2nF4oi/(T2)) = (0,12,5,11,1,7) n A (37) 
A^* = GoF~^(A'^_i) = (1,4m- l,5,4m + 3) n A 3 < m < M. 

The rightmost equahties mean that the polygon on the left can be assembled with an 
appropriate choice of the components listed on the right, in the given order. The issue 
of uniqueness of a polygon with such a labelling is unimportant, as these polygons are 
continuous images of well-defined objects. The boundaries of the sets A^* are specified 
by the following non-empty intersections 

5, 4m - 1 n A^* 7^ m > 1 (38) 

and no other. In particular. A""* is open. The sets (f37l) are indeed quadrilaterals, apart 
from A2'^*, which is a hexagon. 

All atoms are tangent to 1, and we let Po{m) be the left end-point of the intersection of 
A^* with 1. Then we let Pi{m — 1) be the vertex of A^* which is opposite to PQ{m) (the 
term 'opposite' refers to the fact that all atoms have an even number of sides). We find 
that -Pi(O) and i-'i(l) are, respectively, the lower end-point of the intersection of A°^* and 
A2"* with 0; likewise, for m > 2, the point Pi{m) is the left end-point of the intersection 
of A^* with 5^. With these prescription we can characterize the sets AJ^* in terms of the 
vertices Po,i) as follows 

KZ' = [Po{m),PQ{m-l),Pi{m-l),Pi{m)] m > 2. (39) 

A complete description for the cases m = 1, 2 it not needed at this stage, and will be given 
in section 14. 1[ 

We can now state the main result of this section. 

Theorem 11 The regular atoms A^* of the map L°^^ are the polygons ([5^, with the 
boundary specified in \3^) . For each m, the atom A^* is invariant under L°"*, and its 
return symbolic dynamics is 

(1,1,1,(0,1)2-^1,1,1). (40) 

The fixed set of is a segment connecting Po{m) to the opposite vertex Pi{m — 1). In 
addition, there is a (possibly infinite) set of irregular atoms of total area O(A^). 



Proof. From lemma \T0\ it suffices to consider the pre-images of Ti and A', and from 
the rightmost equation in ([32]) . we may replace A' by T2. 

From the proof of lemma [T0\ we see that the part of Ti that maps to A under is 
the ff-symmetrical quadrilateral 

T^nHiT^) = (-3,5^,4,-2) C E. 
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To obtain the first atom we apply Go F to it, which maps it back to A. We obtain 
Af * G G (( -7, 1, 0, -6 ) n G(A)) = ( 8, 7, 1, ) n A, (41) 

as desired. 

We have T2 C H, and so G o F^'*(T2) C A. From the proof of lemma \T0\ we see that 
the part of T2 that maps to A under is the F~'^ o //-symmetrical hexagon 

T2 n o H{T2) = (0, -7\ 5\ -2, 4, -6^) c (0, -2, 4, -6^) = ^-^(a;) (42) 

where A'^ is the first atom of the 2-sector map (lemma [9]) . The atom A2'^* of L°"* is the 
image under G o F~^ of the symmetrised hexagon in (I42p . We obtain 

A°"t = ( 0, 12, 5, 11, 1, 7) n A. (43) 

The rest of T2, namely the triangle 

F-2(A' \ a;) = (0,-62,4) CE 

is covered, apart from a set of measure 0{X^), by the second pre-image of the atoms 
A^, m = 2, . . . , M — 1 of the 2-sector map, from lemma [9l The above inclusion, together 
with lemmas [T] and [3 show that F~^ is continuous on A' \ A'^. Then, from lemma O we 
have 

GoF-^{A'^_^) = (1,4m- 1,5, 4m + 3) n A 

as desired 

Verifying that the polygons defined in (j37p have the boundaries as given in (|23p is a 
simple exercise in orientation, recalling the conventions of table 1. 

The code of each point in A begins with the string (1,1,1,0), associated with the map 
^out _^ (lemma [1]). One verifies directly that the points in A"*^* return to A with the 
code (1, 1, 1). For m > 1, the m-th atom A^* follows the itinerary: 

2 iterations to reach A', with code (1,0) (lemma [7|); 

4(m — 1) iterations to reach H{A'), with code (1,0)2™^^ (lemma[9]); 

4 iterations to reach H', with code (1, 0, 1, 1) (lemma [7]); 

1 iteration to reach A, with code (1). 

Putting everything together, we obtain the code (j40p . In particular, the transit time of 
the m-th atom is odd, and equal to 4m + 3. From lemma HI we then have 

Fix L^"* = (Fix H) n A^"* , (44) 

and it remains to show that all these intersections are non-empty. 

We identify the fixed set of an involution within each regular atom. The atom A°"* 
is the image under G o F~^ of the set Ti n H(Ti), which intersects FixH. Hence A°"* 
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intersects GoF ^(Fixi/) = F^{FixH), in accordance with (dH). Likewise, the atom Ag"* 
is the image under G o F""^ of the set T2 n F'^ o H(T2) which intersects F-^{FmH). 
Hence Aj"* intersects G o F^^{FixH) = F^(FixH), as desired. 

The symmetry properties of the atoms A^* for m > 3 are proved with an argument 
analogous to that used in the proof of theorem \6\ Specificahy, in place of ([26]) and ([27|) , 
we have the symmetric points 

j'^ n H{^'J G Fix n S' m = 0, . . . , L(M - 1)/2J (45) 

as well as 

F^™-2(7^)ni?oF^'"-2(7^) G FixF m = l,...,[(M-l)/2j. (46) 

These symmetry lines intersect images of atoms, and by application of G o -F~*, for a 
suitable t, one ensures that the intersection in (j44p is non empty also for m > 2. We omit 
the details, but we note that, by construction, these sets intersect one vertex of the atom, 
and hence also the opposite vertex, due to symmetry. 

All regular atoms are maximal, being bound by images of the discontinuity line. Then, 
from lemmas [2] and El they are invariant under the map L°"*. D 



4 Quantitative results 

In this section we supplement the geometrical approach to return-map dynamics developed 
in section [3l with a predominantly algebraic one. To illustrate the technique, we will prove 
a version of theorem 111! As we shall see, the proof is conceptually simple, and the details 
relatively easy to carry out, provided that one is willing to rely heavily on computer 
assistance to manipulate complicated expressions. For the sake of brevity, the polygons 
treated in this section do not include their boundaries. 



4.1 Calculation of atom vertices 

The elementary ingredients of our algebraic calculations are the regular atoms of L™ and 
j^out expressed in terms of their vertices. To obtain explicit expressions for the latter as 
rational function in A, we rely on formulae ()22[) and ()37p . which express the polygonal 
atoms in terms of bounding lines. Our strategy is to express the latter in terms of the 
cartesian coordinates of their endpoints, then calculate the intersection of neighbouring 
lines to obtain the corresponding polygonal vertex. All quantities are restricted to ri, that 
is, we don't make use of the periodicity of the torus. 

Before proceeding to the calculation of vertices, let us introduce the convenient notation 
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where i G {0, l}'^ is an (arbitrary) codeword, and 



2 ^0 I X \ ^ I X \ I I 



where the matrix C was defined in ([1]) . From ^ we see that our piecewise isometry F acts 
as F(o) (generahzed rotation about the origin) on the atom Qq and as (generahzed 
rotation about the fixed point (^^a) 2^)) ^i- crucial to distinguish between the 
expressions -f'(ton---n._i); an isometry of the plane, which acts as a generalized rotation 
about some fixed point, and F^, a piecewise isometry of the unit square, acting differently 
on different domains. 

Let us calculate the vertices Pq (m) and Pi (m) of the regular atoms of , see equation 
(j39l) . where, based on the geometric analysis of section [3l 



with 

Cm = -^(I3(0,l){2m-l)l2)([(l,0), (1, 1)]). 

We define 

a = 27r/) = cos-nA/2) = | - q = sin-^(A/2). (47) 

From ([1]) we see that a is the angle of rotation of the matrix C, while = 2tt{^ — p). 
Hence 6 is the angular departure from the A = rotation, that is, 40 is the (positive) 
angle between successive regular components of the sector map on S, in the Q- metric. 



To evaluate explicitly for general m, we regard the isometry P(i3(o,i)(2™-i)i2) as a ro- 
tation by 3a about zi = (2^, jrx)} then a rotation by (4m— 2)a about Z2 = 1^1?) ' 
and finally a rotation by 2a about zi. A multiple-angle rotation corresponds to a power of 
the matrix C and can be executed in a single step, thanks to the formula (a consequence 
of the Jordan decomposition of C) 

N 1 A sin(A;a) / A — 2 \ 
C' = -osika)^^ 1 j+7CT(2 -A J- 

Our geometric constructions of the atom boundaries typically involve powers of C^, in 
which case ()48p can be conveniently written in terms of the 0(A) angle 9 defined in ()47p : 

^41 1 0\ sm{4W) / A -2 \ 

c'' = -osim[, 1 J -71=1(2 -A )■ (^9) 

Once the segments have been determined, the indicated intersections are calculated by 
elementary algebra (a tedious exercise, best done with computer assistance [12]). With M 
as given in (j29]l . we have 

1 ' - ^2 , \4 I ^1 o\2 , \4\ tan^ 



Pi(m) = - 1 + 2A-A^ + A* + (1-3A^ + A*) 



2 V tan(2m - 1) 

A^ A , . n , tan 9 \ , ^ 

l + — + -(X -2) — , 2<m<M. 51) 

2 2^ ^tan(2m - l)6'y ' " " ^ ^ 
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We recall that for m = 1,2, the symmetry axis of A^* intersects instead of 5-^, 
and Pi{m) was defined to be the right-hand endpoint of the symmetry axis (see remarks 
preceding equation (j39p . and theorem llip . The atom Aj"^* is a quadrilateral (see equation 
[37|) . with vertices 

Ar = [^'o(l),(l,l),^i(0),P'], (52) 

where 

The atom A2'^* is the hexagon 

A°"* = [Po(2),Po(l),P",^i(l),P'", A(2)], (53) 

where 

/ 2-2A-4A2 + A3 + A4\ / i _ a - SA^ + + 

■^Ul) = 1'^^; nx2 , ^ ' ^ = 1 



2 - 9A2 + 6A4 - A6 y ' \^ ' 1 _ 6A2 + SA^ - A^ 

1 - 6A2 - 3A^ + llA^ + 4A^ - 12A6 - A^ + 6A^ - A^° 
1 - 6A2 + 5A4 - A6 ' 
1 - A - 4A2 + 3A3 + 4A^ - TA^ - A^ + 5A^ - A^ 
1 - 6A2 + 5A4 - A6 



To calculate the vertices Qo{n) and Qi{n.) of the regular atoms of L™, we apply to 
the line 70 = 1 a 4n-fold C-rotation about the centre of followed by a G-reflection and, 
finally, intersection with the lines 70 and 71 = 5^, respectively. 

The results are [12] 

= [ -A.„(A) + (4-A2)c.(A) 'V ^''^ 
^ , , ^ a(A) + 6(A)^„(A) + c(A)c„(A) f{X) + g{X)sn{X) + hiX)cn{X) \ .... 



d{X)sn{X) + e{X)cn{X) ' d(A)s„(A) + e(A)c„(A) 



with 



a(A) = 2 - A - 7A2 + 3A3 + SA^ - A^ - A^ 6(A) = -1 - 2A - A2 + 2X^ + A^ 

c(A) = 2 + A - 6A2 - 3A3 + 2X^ + A^ d{X) = -A(5 - 5A2 + A^) 

e(A) = (4 - A2) (1 - 3A2 + A^) /(A) = -A(4 - 2A - 4A2 + A^ + A^) 

g{\) = A(-3 - A + 2A2 + A3) h{X) = 4 - 7A2 - 3A3 + 2A^ + A^ 



and 



Cfc(A) = cos{4k9) = T^k 
Sk{X) = V4 - A2 sin(4A;0) = (^^) U,k-i 
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Here and Uk are the Chebyshev polynomials of the first and second kinds, respectively 
[19, p 1032]. 

The use of the rotation angle 4n0 in (j49p . as well as in the formulae for Qo{n) and 
Qi{n), allows for a natural interpolation of these quantities in the range < n < — 1. 
For continuously varying n, the matrix C^"" defined by (j49p (not (j48p ) represents a family of 
C-rotations whose angle increases monotonically from to 4:N9 = tt/2 — 0(A). Moreover, 
as n increases from to — 1, Qoin) and Qi{n) smoothly trace out, in order-preserving 
fashion, the upper and lower boundaries of the regular part of A. 

Note that formulae (j54p and ()55p explicitly specify the invertible transformations be- 
tween parameter n and points X = {x,y) on the lines 70 and 71, respectively. This gives 
precise meaning to the functional notation n{X),X £ S, which we will employ below. In 
section [5] we will exploit this parametrization to study the relative positions on the two 
boundary lines of the points Po^m) and Qo{n) (resp. Po{in) and Qo{n)). 

4.2 A preliminary lemma 

The following lemma plays an important role in the theorem which follows. 

Lemma 12 The map is an isometry on the triangle [(0, 1), (1, 1), (A^, 1 — A)], with 
itinerary 

[(0, 1), (1, 1), (A2, 1 - A)] ^ [(0, 0), (0, 1), (A, A^)] ^ [(0, 0), (1, A), (1, 0)] 
^[(0,0),(0,1),(A,1)]^[(0,0),(A2,A),(1,0)]. 

The triangles F^{[{0, 1), (1, 1), (A^, 1 — A)]), A; = I, ... ,4, do not intersect the sector S. 

Proof. The proof is by direct mapping of the vertices, accompanied by straightforward 
verification that the interior each triangle is contained in a single atom of F, the latter 
being correctly identified in the statement of the lemma. The empty intersection of the 
four image triangles with S is trivially verified. 

4.3 Algebraic alternative to theorem 1111 

We now state an alternative to theoremllllin which the polygons A^* are defined explicitly 
in terms of their respective vertices, as in section 14.11 That these are symmetric atoms of 
2^out -g p];.Qvgfj a,n inductive calculation of their return paths. 

Theorem 13 For 1 < m < M , we have 
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Moreover, 

[PQ{m),Pi{m - 1)] = FixL°"* 



Proof. We begin by defining some auxiliary polygons related to A^*: 
Ar = U^L Ar U [Po(M), Pi(M), {g, 1)], 

^OUt d£f i?|,^3^g^-^^2m-l) o G(A5^*). 

We want to prove that G(Am) is mapped isometrically onto A^ by F^^^^ . 

The proof is by induction on m. For m = 1, the isometric property follows from lemma 
[12] by showing that ^(GA^"*) is contained in [(0, 1), (1, 1), (A^, 1 - A)]. To verify this, we 
note that 

F{GkT') C F{[{g'\ 1), (1, 1), (1, g)]) = [(A/, g"), (A, 1), (1 + A - 5, !)]• 

It is a simple exercise [12] (sec. E.5.1)to show that for < A < 1, all three vertices of the 
image triangle are within or on the boundary of [(0, 1), (1, 1), (A^, 1 - A)]. 

We now assume that the induction hypothesis F^m+i ^ (^(yVm) = A^ holds for a given 
m < M. Note that G{Am) is partitioned by the line segment G{[Po{m), Pi{m)]) (m > 1) 
or G{[Po{l),P']) {m = 1) into adjacent polygons G(A^*) and G{Am+i)- We denote the 
images of these polygons A+ and A~, respectively. One checks [12] (sees. E.5.2, E.5.3) 
that the boundary separating the two is a subset of [(0,0), (1, A)], so that 

A- c no, A+ c ([(0, 0), (A2, A), (1, 0)] \ no) C (1^1 \ S). 
To complete the proof of 



-i4m+l 



G(A„+i) = A„+i, (56) 



we need to show that maps A^ isometrically onto Am+i, with itinerary (0,1,0,1). Using 
lemma [T2] and A~ C ilo; we have 

F2(A,-) = F(oi)(Ajcf)o, (57) 
where the inclusion is verified explicitly [12] (sec. E.5.4). Again using lemma [T2\ we have 

^^A") =F(oioi)(A„). 
The right-hand side is identical to A^+i by construction. 

Having established (j56p . we need only show that F^ maps A+ isometrically onto A^* 
with itinerary (1,1). For the isometry, we verify the inclusion [12] (sec. E.5.5) 

Fi(A+)c(l^i\S). (58) 
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The rest is an explicit calculation [12] (sec. E.6) of the image points of the vertices, using 
elementary algebra, In particular, we find that the map LJ^* = ^F^m+s ^ q ^^^^g ^j^g 
vertices of the polygons A^* as follows (notation of section HTT]) : 

Po{m)^PQ{m), Pi(m - 1) Pi(m - 1), 1 < m < M, 

Po{m - 1) ^ Pi{m), 2<m< M, 

(1, 1) ^ P', m = 1, P" ^ P'", m = 2. 

The final identification of _p4m+3 o G as the first-return map , with regular atoms 
A^* requires verification that all polygons oG{A^^), 1 < k < 4m + 3, are disjoint from 
E. We have taken care to establish this property at each stage of the proof (specifically, 
in lemma [12] and equations ()57p and (I58p ). 

Recalling that the map acts linearly on A^*, the vertex mappings are sufficient 
to establish that the atom is L°"*-symmetric with fixed segment [Po{m), Pi{m — 1)]. D 



5 Regular atoms and fixed points of L 
5.1 Overview 

In section [3] we decomposed the return map L into a product of involutions L™ and 
L°"*. According to theorem [6l is endowed with a sequence of regular atoms A™, n = 
1, . . . , N — 1, each of which is an L™-symmetric polygon. Theorem [TT] establishes an 
analogous sequence of regular atoms, AJ^*, m = 1, . . . , M for We now define regular 
atoms of L to be the nontrivial intersections of the regular atoms of the two involutions: 

A™,„ = AT n Ajf . (59) 

For some, but not necessarily all, of these atoms, the symmetry lines of and L°'^^ in 
A™ and AJ^*, respectively, intersect within Am,n, at a symmetric return- map fixed point 
Z{m,n). Such fixed points will be called regular, and will be the main focus of this and 
succeeding sections. They correspond to a 2-parameter family of periodic orbits of F, with 
periods 

t{m, n) = 4(m + n) — 1 (60) 

and symbolic codes 

i^"-'") = (l4"-i,(0,l)2™-i,l,l), (61) 

the latter obtained by concatenating the codes of the involutions L™'°"* given in theorems 
Elandllll 

Before examining in detail the geometric properties of the regular atoms of L, it is 
useful to introduce a broad classification according to the size of the index m relative to 
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Figure 7: Top: Superposition of the atoms of the domain A into atoms. The atoms of L™ are in 
red; those of in black. Bottom: Magnification of the region near the tip of A, showing the 
crossover phenomenon. The integers denote values of to, and crossover occurs for m — Q. 

its maximum value M = M{X), given by equation (j29p . A typical example is displayed 
in figures [7] and m where A has been assigned the value Figure [7] shows the regular 
atoms of (lighter boundaries) and (darker boundaries) in all of A (note that the 
aspect ratio has purposely been distorted to make the structural relationships visible). 
The bottom figure \7\ zooms in on the left-hand end of A, while figure [8] zooms in on the 
atoms with m = 3, 30 < n < 40, showing explicitly FixL™ and FixL°"* and the regular 
fixed points of the return-map where these symmetry lines intersect. 

We observe some interesting features of the regular atoms and their symmetry lines. 
The orientations of the 'vertical' boundary lines change substantially as one traverses A. 
The A™ boundaries are vertical at the right-hand edge but suffer significant counterclock- 
wise rotations as one moves toward the left-hand vertex. The A°"* boundaries, on the 
other hand, rotate in the opposite sense. Somewhere in the middle, there is a crossover 
region where the boundaries of the two sets of atoms are approximately aligned (this cor- 
responds to m = 6 in figure [7]). To the right of the crossover region, each intersects 
a number of successive A™ atoms, and each of the intersections Am,n (the regular atoms 
of L) contains a unique regular fixed point Z{m,n). To the left of the crossover region, 
on the other hand, each A™ intersects a number of successive A°"* atoms, and now there 
are regular atoms A^.n which do not contain a regular fixed point. This can be seen quite 
clearly in the case m = 8, which has a triangular intersection region in its northeast corner 
which is not intersected by the diagonal Fix Lg"* . 

In appendix A. 7 we give a precise definition of the crossover point m = m*,n = n* 
and obtain a perturbative formula for m* and n* as functions of A. In particular, we will 
find 

m = —= 1-F-H , n = — 1= - T H • 

V2A V 3 y ' 4A V2A 4 

We note that in the limit A ^ O"*", the crossover value of m tends to infinity. Thus, if 
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Figure 8: The fixed points of tlie return map L are tlie intersection of tlie fixed lines of the 
involutions 2^'"^°"*^ which are segments connecting opposite vertices of the atoms A'" (vertices 
QoW, Qo(n - 1), - 1), Qi{n)), and A^"* (vertices Po(m), Po(m - 1), Pi(m - 1), Pi(m)). The 
example focuses on m = 3, n = 34 for A = 

we are interested in the asymptotic behaviour of regular fixed points and their stability 
ellipses for fixed positive integer m, we can limit our attention to those regular atoms 
which are well to the right of the crossover, and avoid the complexities of the crossover 
phenomenon. This is achieved by introducing a cut-off condition 

m < X-""", where < z^o < ^- (62) 

We will prove that the regular atoms satisfying this condition have a number of nice 
properties, which allow us to obtain explicit formulae for their total area, whereas those 
with m exceeding the cut-off have a total area tending to zero as Aa"*^", and can therefore 
be neglected. 

In the remainder of this section, we study in detail the properties of the regular atoms 
of L. Our starting point is the list of explicit formulae for the atomic vertices given in 
section |4?TJ The vertex coordinates are rational functions of A; as A ^ 0, we will consider 
Taylor expansions in A, with suitable remainder. Our most complete description of the 
atoms Am,n will be restricted to those satisfying cut-off conditions A < Aq, m < A~'^". Our 
choice of cut-off parameters is (cf. (j62|) ) 

Ao = 10"^ and 1^0 = ^- (63) 

Even though a specific choice is required to obtain explicit estimates, our qualitative 
results are insensitive to these special values, and the detailed calculations included in the 
Electronic Supplement [12] can readily be adapted to other cut-off values. For the rest of 
this section, we shall assume that Aq and fo are given by (j63p . 

The organization of the regular atoms of L™'°"* and their intersections Am,n is deter- 
mined by the order and spacing of the vertices along the lines 1 and 5^. Often it will 
be convenient to regard n as a continuous variable on the interval [0, A^], as discussed in 
section 14. H and to express our results in terms of auxiliary variables 

TT 

?7 = — — 2n9, T = tan(?7). (64) 
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Note that as n increases from to N, the quantity r decreases monotonically from 1 to 
tan(7r/4 - INO) = 0(A). 

In section 14.1^ we introduced the notation n{X) to designate the real parameter n 
assigned to a point X on one of the hues 1 or 5^. In the same spirit, we will employ t{X) 
to denote the corresponding r- variable, namely tan(7r/4 — 2n{X)9). 



5.2 Expansions for the atom vertices 



We begin our structural analysis by determining the positions of the vertices Po,l("^), Qo,i{n), 
correct to next-to-leading order in powers of A. The following lemma establishes, among 
other things, that the horizontal positions and widths of the atoms AJ^* become fixed as A 
tends to zero, while the vertical dimension shrinks proportional to A. For the atoms A™, 
on the other hand, both dimensions shrink proportional to A. In addition, we see that the 
j;-coordinates of Po{m — 1) and Pi{m) coincide up to first order in A. These properties 
are apparent in the figures. The proof of the lemma is given in appendices A. 2 and A. 3. 



Lemma 14 For A < Aq, and 2 < m < X , with Ao,fo given in \63\). the cartesian 
coordinates of the vertices of the regular atoms A^* and A™ are given by 

m + 1 



Po{m) 
Pi{m) 

Qo{n) 
Qi{n) 



2m + 1 

m 
2m -1 

1 + 



+ A 



2m + 1 



6 



+ A-^fi^^^±^K^^A^ + n.A^), 



2m - 1 



6 



1 



2m - 1 



^^m(m + 2)^3^^^^^. 



1(1 + r) + ^(3 - 2r - r2) + ^(2 - 5t + 2t' + r^) + r^.A^, 1 
1 A A^ 

-(1 + r) + -(7 - 2r - Sr^) + — (-6 - 29r + lOr^ + 25t^) + rg.A^ 
Z o oZ 

A3 



(65) 



(66) 



\2 \3 \ 

1 - tA + — (1 + 2t + 5t2) - —(2 + llr + lOr^ + 25t3) + r^yX^ , 
4 16 / 

where r is defined in (^Jj, and the remainders r are rational functions of X^^"^ that admit 
the following bounds 



r2x G 



1697 1507 



IQS ' 105 
2051 4499 



5000 10000 

In addition, we have 



Tlx G 
r2.x G 



8523 1243 



106 ' 105 
23 389' 
20' 250 



riy G 

G 



5449 1203 



106 ' 50000 
519 539^ 



Po{l) = (l + X-lX^ + roixX\ 1 
Pi(0) = (l,l-iA + |A2-A3 + rioj,A4), 
Pi(l) = (l, 1 - A + |A2 - 4A3 + niyX^ 



5000 100 



463 3703 



roix G [ 12500' 100000 J 
r a [1249 313] 

' 1% ^ L 500 ' 125J 

r a [ 8739 8763 ] 
'lly ^ LlOOO' 1000 J ■ 
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Corollary 15 For A < Aq and 2 < m < A '^'^ , we have 

Poim)^ < Pi{m)^ < Po{m - 1)^ < Pi(m - 1)^, 

with 

Pi{m), - Po{m), = 0(1), Pi(m), - Po(m - 1), = 0{X^). 

For m = 1, 

Po(l)<A(0), Pi(0)-Po(l) = 0(l). 

In addition to the information about the horizontal ordering of cartesian coordinates, 
we will need analogous results with respect to n, considered as a continuous variable. The 
proof is found in appendix A. 5. 

Lemma 16 For A < Aq and 2 < m < X^'^° , we have 

n{Po{m)) > n{Pi{m)) > n{Po{m - 1)) > n{Pi{m - 1)), (67) 

and the quantities X{n{Po{m))—n{Pi{m))) and X^^{n{Pi{m))—n{PQ{m—l))) have positive 
lower and upper bounds which are independent of X. In the case m = 1, we have 

n(Po(l)) >ri(Pi(0)) =0, 

and the quantity A(n(Po(l)) — "'(-fi(O))) has positive lower and upper bounds which are 
independent of X. 

5.3 Structural theorem 

We are now in a position to establish the most important structural properties of the 
regular atoms of the return-map L. 

Theorem 17 If the inequalities 

Poiin)x < Qoin - l)x and Pi(m - 1)^; > Qi(n)a; (68) 

are satisfied for some integers m, n in the range 1 <m < M , 1 < n < N , then the regular 
atom Am,n = ■'^Jfn A™* has positive area, and contains a regular fixed point Z{m,n). The 
stability ellipse £rn,n surrounding Z{m,n) is tangent to (at least) three sides of Am,n, one 
of which is a subset of either 1 or 5^ . 

Given < A < Aq and 1 < m < A"''", then there is a non-empty set of values of n for 
which the inequalities Ii68\) hold, given by 

n-(m) < n < n^{m), n_(m) = [n(Pi(m — 1))] , n^{m) = \n{P(){m))~\. 

As X ^ 0, such a set becomes infinite. 
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Proof. From equations and ([39]) we see that the points Po{m) and Qoin — 1) 
belong to the segment 1, which is parallel to the x-axis. Suppose that m > 2. Then the 
points Pi{m — 1) and Qi{n) belong to 5^, which is quasi-parallel to the x-axis. Then the 
inequalities (168]) imply a corresponding ordering of the Ps and Qs along the respective 
segments. It follows that the segments [Po{m), Pi{in— 1)] and [Qo{n— l),Qi{n)] intersect 
transversally at a point Z (171,71), which lies in the interior of both atoms, hence in the 
interior of their intersection. This establishes the existence of a regular atom Arn,n of 
positive measure. From theorems iGlandfTDwe recognize [Po{m), Pi(m— 1)] as the symmetry 
axis for L^^, and [Qo{n — 1), Qi{n)] as the symmetry axis for L™, so that the intersection 
point Z{m, n) is indeed a regular fixed point of L. If m = 1, 2, the symmetry axis of L°^* 
connects -Po("i) to 0; the second inequality in (f68]) is automatically satisfied, and the first 
inequality in (|68p suffices to establish transversal intersection of symmetry lines. 

Because the atoms AJ^* and are maximal, so is A^ jj. Since the fixed point Zijn, n) 
lies in the interior of Am,n, the associated ellipse Sm,n is necessarily tangent to the boundary 
of Am,n- In particular, Z{m,n) is stable. 

We must now establish the number of points of tangency. Let T be a point at which 
Am,n is tangent to £m,n- To make our argument independent from the details of the 
boundary, we use the expression Lm,n"*(r) to denote the limit of L^m°n^{z) as z approaches 
T from the interior of the atom. Now, any L-invariant ellipse in the interior of £m,n is also 
invariant under the individual involution and LJn*- To see this, note that, generically 
(irrational rotation number), such an ellipse contains a symmetric dense orbit, which is 
mapped into itself by one involution, hence by the other. By continuity, this result then 
extends to the non-generic case (rational rotations). Because of this property, if T is 
tangent to £, so are L™{T) and L^^{T). 

Let us assume that {m,n) 7^ (1,1) or (2,1), which are special cases to be dealt with 
separately. We distinguish three possibilities, depending on whether A^,™ is tangent to 
5, to 1, or to botlH (see figured]). If A^.n is tangent only to 5, then T —4n + 4, 
because L™(— 4n -|- 4) = 1, which lies outside the closure of the atom. The two involutions 
exchange the remaining three boundaries as follows 

— 4n <— ^ 5 4m — 1 

and since the images of T lie on the boundary of £, we get one tangency point on each 
of the above three sides. In particular, if Arn,n is a triangle, then £rn,n is tangent to all 
its sides. In the case {m,n) = (2, 1), the atom A2,i has the additional involutory relation 
12 = U^^{0). Since T cannot be tangent to 0, then it cannot be tangent to 12. But then 
A2,i is necessarily a pentagon (i.e., 5 is a side of the atom), for otherwise the remaining 
two sides would map outside the atom, and we would have no tangency, contradicting 
maximality. So we still get three points on 5, 7 and 12. 

An analogous argument applies to the case in which Am,n is tangent to 1 but not to 
5. Now the ellipse £ cannot be tangent to — 4n, and the involutions act on the remaining 

^To lighten up the notation, we'U omit aU superscripts identifying regular components. 
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sides as follows 



r in T out 

-4n + 4 ^ 1 ^ 4m + 3 



giving again one tangency point on each of the above three sides. In the special cases 
(m, n) = (1,1), we have the additional involutory relation 8 = L°"*'(0). Again, Ai^i is 
necessarily a pentagon (lest we would have no tangency). This time the ellipse is 
tangent to 0, and so we obtain a fourth point of tangency to 8. 

It remains to consider the case in which A^.n is tangent to both 1 and 5. We cannot 
have the quadrilateral atom Am,n 7^ A™, because its boundaries 1 and 5 are mapped 
outside the closure of the atom by L^*' ^"^^ mapped to the other two boundaries by LJf. 
But this would imply the absence of tangency points, contrary to the maximality of the 
atom. For a similar reason Am,n 7^ ^m*- Thus Am,n is either a pentagon or a hexagon. 
In the former case, symmetry still forces the same three tangency points as above, one 
of which is either 1 or 5. However, £ cannot be tangent to both, since one of these 
sides is mapped outside the atom by L^*- the hexagonal case there are two possible 
configurations with three tangency points, of which at least one must be realized. Each 
configuration must be realized in some parametric interval, and by changing parameter the 
two configurations exchange role. Then there must exist a parameter value corresponding 
to six tangency points, the most symmetric configuration. 

Let us now assume that the parameters A, m, and n belong to the given range. Suppose 
-Po("^)x > Q^in — Then T{Po{m)) > T{QQ{n — 1)). For m = 1, the fact that Qi{n — 1) 
lies on 5 implies that the right-hand edge [Qo{n — l),Qi{n — 1)] of AJf lies to the left of 
the interior of A°"* and there is no nontrivial intersection. Assume therefore m > 2. Since 
T{Qo{n- — 1)) = T{Qi{n — 1)), and, using ([67j) and monotonicity, T(Pi(m)) > T{Po{m)), 
we conclude that T(Pi(m)) > T{Qi{n — 1)). Thus [QQ{n — 1), Qi{n — 1)] lies to the left of 
the interior of A^* and there is no nontrivial intersection. 

If m = 1, the second inequality in (j68p is automatically satisfied, since Pi{Q)x = 1- 
Assume therefore m > 2, and suppose Pi{m — l)^ < Qi{n)x- Then r(Pi(m- — 1)) > 
r((5i(n)). But T{Qi{n)) = t((5o("-)) and, using ([67|) and monotonicity, r(Pi(m — 1)) > 
T{Po{m — 1)). Thus r(Po("^ — 1)) < T^Qoi^))- We conclude that the left-hand edge 
[Qo{n'),Qi{n)] of A™ lies to the right of the interior of A^* and there is no nontrivial 
intersection. 

Defining n+(m) = \n{PQ{m))~\, we know that n^{m) — 1 < n{PQ{m)) and hence, 
by monotonicity, the first inequality of (j68l) holds, not only for n = n^(m), but for all 
n < n+{m). Similarly, if n-{m) = \n{Pi{m — 1)], we know that n-{m) < n{Pi{m — 1)) 
and hence, by monotonicity, the second inequality of (|68p holds, not only for n = n-.{m), 
but for all n > n-{m). We conclude that both inequalities hold for n-{m) < n < n+(m). 

From lemma [T6| we know that X{n{Po{m)) — n[Pi[m — 1)) is bounded below by a 
positive number. Hence n^{m) — n-{m), which by definition exceeds n{Po{m) — n{Pi{m — 
1)) — 1, tends to infinity in the limit of vanishing A. D 
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5.4 Expansion of Z{m, n) 



The point Z{m,n) lies at the intersection of the hne segments [PQ{m), Pi{m — 1)] and 
[Qo{n — l),Qi{n)]. Using hnear algebra, we have 



Z{m, n) 



1 



/ -Zx 






X 


w \ 


V -^y 


Wy ) 


(: 


X 


^ J 



(69) 



where 

u = Po("i), v = Qo{n-l), w = Pi{m - 1) - Po{m) , z = Qi{n) - Qo{n - 1). 
and, for any pair of 2-vectors a and b, the cross product is defined as 

a X b = Uxby — ttybx- 

Inserting the estimates of lemma [TJ] and corollary and once again applying the rules 
of interval arithmetic, we now obtain similar estimates for Z{m,n). 



Lemma 18 For < A < Aq and I <m< X''"" , 



Z{m, n) 



1 



A 



1 + r) + - (7 + (1 - 4m)r2) 



\2 

+— (4 - (5 + 16m)r + (1 - 8m + IQui^^) + rzx^^'^. 

A A^ 
! + -(!- (2m + l)r) + — (l + 2m + (8m^ + 2m - 1)t^) 

\3 

H (-20 + 48m - 16m^ + (23 + 6m - mm?)T - (3 - 18m + 96m^)r^) + rz 

192 



2000 ' 5000 



where t = T{QQ(n)) and, for m > 3. 

1033 2671 

rzx G 

For m = 1,2, we have 
Z(l,n) 

with 



rzy G 



5563 7373 



1000 1000 



l{l + r) + ^{7- 3r2) + ^(4 - 21r + 9t') + rzixX\ 
1 + -(1 - 3r) + ^(1 + 3r2) + _(4 - 17r - 27t^) + rziyX'j , 



rzix G 



1011 334 



and 



Z(2,n) 



100 25 



7A, 



rziy e 



4211 437 



100 10 



;i + r) + ^(1 - t2) + ^(4 - 37r + 49r3) + rz2xX\ 



\ c\2 \3 \ 

1 + -(1 - 5t) + ^(1 + 7r') + -(4 - 95r - 245r3) + rz2yX'j 
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with 



rz2x e 



506 331" 



rz2y G 



9517 

loo 



,150 



The proof is given in appendix A. 6. 



5.5 Disk areas 



In this section we complete the proof of theorem B. The regular fixed point Z^m, n) of L is 
a periodic point of the map F, with period t{m, n) = 4(m + n) — 1, see equation (j60|) . Let 
-A-m,n be the area of the cell (stability ellipse) of Z{m, n); furthermore, for fixed m, let Am 
be the total area of the ellipses associated with orbits of Z{m, n), n = n^{m), . . . , n-|-(m). 
Thus 

n_|. (m) 

Am= ^ Am,nt{m,n) t{m,n) = 4:{m + n) - I. 

n=n- im) 

Finally, we denote by A^'^^ the total area associated with the regular orbits: 

M 



771=1 



Our goal in this section is to prove that the A — > 0"*" limits of A^a and A^'^^ exist, and 
to derive explicit expressions for the limiting areas. 

From theorem 1171 we know that the ellipse associated with Z{m,n) is tangent to one 
of the lines 1 or 5^, and so the area of the ellipse is just vr times the square of the Q- metric 
distance from Z{m, n) to the tangent point on the closer line. Using the inner product Q 
defined in equation ([1]), the distance between a point Z and a line segment [P, Q] is given 

by 

^mm^^^y U = Z-P,V = Q-P. 

For the area of disks tangent to the line 1, we apply this formula with Z = Z{m, n),P = 
(0,1), Q = (1, 1) to obtain the expression 

A(Z = ^{l-^yZ{m,n)y-lf. (71) 

For the disks tangent to 5^, on the other hand, we use P = {g, 1), Q = (1, g") (cf. equations 
(fTTI) and (mi)), to obtain (writing Z = {Z^,Zy)) 

A^l = ^ (l - ^) (1 - + (1 - 2Z^)\ - (1 - 3Z,)A2 - (1 - Z,)A3 - ZyX^f . (72) 
The above formulae give us a straightforward way of calculating 
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for given m, n, A. Supplemented by Taylor expansion with rigorous bounds on the remain- 
ders, the same formulae also provide the machinery for our asymptotic calculations in the 
limit A ^ 0+. 



5.6 Bounds for the disk areas 

To obtain formulae for the disk areas (m,n) and A^^\m,n), we insert Z{m,n) from 
lemma [18] into (j7ip and (j72p . respectively, combining the remainder intervals using the 
rules of interval arithmetic [12] (sees. E.10,E.10A). We state the results as a lemma. 



Lemma 19 For < A < Aq and 1 < m < A , 



A 



m,n 
(1) 



^(l-(2m + l)T)2A2 + rAoA5/2^ 
^{l-{2m-3)TfX^ + rAi\'^^ 



(73) 
(74) 



where r = T{QQ{n)) and, for m >3, 



6347 183 



100000 6250 



7091 2857 



100000 ' 50000 



For m = 1, only A^^\l,n) is relevant, with 

221 1179 



rAO G 



6250' 100000 



and, for m = 2, 

393 787 



rAO G 



2000 20000 



rAi G 



7859 7859 



1000000 1000000 



rAi - rAO G 



859 3761 



50000' 100000 



rAi - rAO G 



63 943 



2000 5000 



For fixed m, the fixed points Z{m, n) are lined up along the line segment [Po(^)) -Pi {rn— 
1)] in order of decreasing n, with Z(m, nj^{m)) just to the right of Po{m) and Z{m, n-{m)) 
just to the left of Pi{m — 1) (see figure [9|). Regarded as a continuous function of n, 

^(0) 

(m, n) increases from zero monotonically with decreasing n while ^4^^^ (m, n) decreases 
to zero monotonically. Thus there exists a value n = n^idim) such that 



^(°)(m,nmid("i)) = A'-^\m,n^id{m)), 



and hence 



mm 



{A'-^\m,n),A'-^\m,n)} 



A^^\m,n) n^{m) <n< nQ{m) — 1 



where 



I A^^\m,n) no{m) < n < n+(m) 
no(m) = [nmid(m)]. 



(75) 
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Ql(«mid. 



Ol(no-l)" 




^i('«-1Tq^(„__1) 



Figure 9: Sketches of the atoms A'f^, AJf^, and A™_. The real number rimid is defined by the 
condition that Z{m,nmid) lies on the 'mid-line' equidistant (in the Q-metric) from the upper and 
lower boundaries. 

To locate r^aid = tan(^ — 2njnid^) (see equation (fM|) ). we equate A^^\m,n) and 
A^^\m,n) in [73|) and ([71|) and solve the resulting quadratic equation in r to get (dis- 
carding the irrelevant root) 



^mid 



1 + ^1 + f (2m - l)Ai/2(r^, _ tao) 
2(2m - 1) 



Next we have 



Lemma 20 The quantity A^^\m,n) — A^^\m,n), considered as a function of t for fixed 
m G [2,M], vanishes at 



where, for m > 3, 
and, for m = 2, 



''"mid 



''"mid e 



Im — 1 

2189 599 



^mid £ 



100000 25000 
1003 601 



25000 5000 



Lemma [20] is easily proved using the inequalities (valid for < x, y, < 1 and —x < z < y) 



l-x< Vl + z< l + y/2. 
An immediate (and useful) consequence of this lemma are the bounds 



Tmid < 



5' 



(76) 
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|(2m-l±2)r^id-l| < 1, (77) 
|AW(m,n)| < ^A^, n„,id<n<n+, (78) 

|^(^)(m,n)| < n_<n<n„,id. (79) 



5.7 Total disk area 

We are now in a position to write an exact formula for the sum of disk areas within a 
single atom A'j^^: 

n+{l) 

A, = Yl A^''Hhn)t{l,n), (80) 

n=l 

no(m) — 1 n_|_(m) 

= ^ A(^)(m,n)t(m,n) + ^ A^'^\m,n)t{m,n), (81) 

?i=n_(m) n=n.o('^) 

where the 'midpoint' no{m) was defined in (j75p . To calculate the leading behaviour for A 
tending to zero, we note that the terms in (j8ip proportional to 4m — 1 are negligible in 
the limit. Specifically, 

3^1 < 3Af^(°Vl,l) < 12iVA2 < — A^/^ (82) 

25 

(4m -1)^^ < 4mA^y4(°)(m,nmid) < ^A^/^ m > 2, (83) 
where we have applied ([78]) as well as 

A < 10-^ m < A-l/^ AiV < ^ < 1. 



The next step toward the derivation of an asymptotic formula is to insert the estimates 
(1751) and (HH) of lemma [ig to obtain 

^1 = ^ E (l-3r)^n +r^AV^ (84) 



n=l 

7io(m) — 1 n+(m) 



= ^ I E (1-(2"i-3)t)2 n+ (1 - (2m + l)r)% 1 + r^A^/^ 

yn=n_(m) n=no(m) y 

where r = r((5o(^)) and 

, , 1 
FA < 7- 
5 

Next, we replace the sum over n by an integral over rj using the Euler-Maclaurin formula: 
Ai = — (inn(l-3T)2 + — (4 + (l-3tan?7+)V(l)) + (™o + rA)A^/^ 
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-4n 



{I - {2m - 3)Tf ndn) + (/ {I - {2m + l)Tf ndn 



no 



+ — ( (1 - {2m - 3) tan r?_)^ n_ + (1 - (2m - 3) tan(77o - 29))'^ (no - i; 
8 

H — — ( (1 — (2rn, + 1) tan r^o)^ no + (1 — (2m + 1) tan r/_(_)^ n^ 



(85) 



+{rEMO + rEMi + ?^a)A^/^, 
where, for 1 < m < A"^/'*, 

TEMO < ^ dn I ^ (n (1 - (2m + 1) tan(f - 2n0)) 

"^''^ f /n""^ (1 - + 1)7")^ + 46'(2m + 1) dnn\l- {2m + l)r| sec^ rj 



< 



8 



< ^(iV + 49N\2m + 1)) < ^ 
and an analogous expression for temi for m > 1. 



As a final simplification, we change the integration variable from n to r/, defined in 
equation (j64p . changing the integration limits to the values ry = cot~^(2m + l), cot~^(2m — 
1), and cot~^(2m — 3). We state the result in the form of a lemma, which refers to the 
cut-off parameter values ([63]) . 



Lemma 21 For < A < Aq, we have 



■Am. 



I{2,l,2)+rX^/\ 

I{m + 1, m, m + 1) + I{m, m — l,m — 1) + rA^^^, 2 < m < A" 



where 



/■cot-i(26-l) 

I{a,b,c) = - / {{2c-l)tan7]-lf{7r/4-7])dr] 

4 Jcot-i(2a-l) 



and |r| < ^. 



Proof. The contribution to r of the remainder terms in (j85p is bounded in magnitude 
by 1/2. Moreover, the endpoint terms in ()85p . for both m = 1 and m > 1, are bounded 
by ^iVA^ and so (using A'"A < 1) the magnitude of their contribution to r is bounded by 
1/50. 

The change of variables produces a factor 

^ = ^ e [1 - ^V3, 1] C [1 - 10-^AV^ 1]. (86) 

Since each of the integrals over r] is bounded by ^ (^)^ < ^, we see that the total contri- 
bution of the variable change to r is bounded in magnitude by 10~^. 
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Finally, we can estimate the shifts in the integration limits using 



-1 



-1 



-1 



?7+ — tan 
??_ — tan 
r/o — tan 
rio + 26 - tan"^ 

Using lemmas [20] and [2 



2m4 
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1 




2m- 
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1 




2m- 
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1 




2m- 


1 



< 
< 
< 
< 



tan-i(r((Po("i)) - taW^ 
tan^^(r(Pi(m — 1)) — tan 
tan^^(rini(j) — tan^-*^ 
tan~^(rini(j) — tan^-*^ 



1 



2m+l 



+ 26, 



-1 



1 



2m-l 



1 



2m-l 



1 



2m-3 
+ 29, 

+ 26. 



26, 



I and the inequality (valid for x > 0, y > 0) 



tan ""^(x + y) — tan ""^(x) = tan 



1 + x{x + y) 



< tan ^(y) < y. 



it is straightforward to calculate bounds for the right-hand expressions. We find that he 
contribution to r of the shifts in the integration limits is bounded by ^. Adding up the 
various contributions to \r\, we get the stated upper bound. □ 



Now we are ready to sum up the disk areas from the orbits of all regular fixed points, 
thus completing the proof of theorem B. 

Theorem 22 In the limit A 0^, the total area A^°^{X) of the islands of the regular 
periodic orbits of F converges to the positive quantity 

oo 

yl'-^g(O) = I{2, 1,2) + ^ (/(m + l,m,m + l) + I{m, m - 1, m - 1)) . 

m=2 



Proof. It suffices to show that 

A'^^iX) = 1(2, 1,2) + ^ (/(m + 1, m, m + 1) + I{m, m - 1, m - 1)) + 0(A^/^). 

m=2 

where I is given by lemma [21] and /i(A) = [A^^/^J . The partial sum 

m(a) 

m=2 

is already in the desired form, since the integrals are A-independent and we have [A~^/'^J 
remainder terms, each bounded by the same constant times A^/^. 

The remaining terms, 

M(A) 
m=fi{X)+l 

are bounded by the (Q-metric) area of the triangle 

[P,Po(MA)),Pi(MA))], P=i ^ ,l] , 
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multiplied by the longest return time, 4(M(A) + N[X)) — 1. Inserting our estimates for 
the atomic vertices, the triangular area is 

5V'1-AV4 (Po(m(A)U - P,) (1 - P,{li(\)),} 

while the upper bound on the period is 0{X~^). Thus the upper bound on the sum of 
terms with m > iJ,{X) is 0(A^/^). □ 



5.8 Numerical evaluation of area formula 



Numerical integration gives the following results: 

Area Am 



m 



0.04394252102495575454 

0.027915747684440153071 

0.0043390573122902285760 

0.0011842200753144612564 

0.00044544206984605774866 

0.00020336198417268636974 

0.00010566669827864850788 

0.0000602301047692083869367 



For asymptotically large m, the integrands can be expanded in powers of ij and integrated 
term by term to yield 



A 

~ 48m4 



+ 



7r(7r — 1) 
24m5 



Summing the numerical integrals over the first 2000 values of m gives 

2000 



m=l 



0.0783220277996. 



This is approximately 36 percent of the area in the square outside the fixed-point disk. 



6 Covering 



In the previous sections we defined the return map L of the domain A, and then we 
identified a prominent sub-domain of A, namely the union of the regular atoms of L. In 
the following theorem we show that the set of regular atoms is an appropriate surface of 
section for the orbits of F oTitside the main island £. In so doing, we will complete the 
proof of theorem A in the introduction. 
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Theorem 23 The images under F of the regular atoms of the return map L cover all of 
0,\£, apart from a set of area O(A^). 



Proof. All metric considerations refer to the Q-metric (j4]). We begin to deal with the 
area of the region complementary to £. The Q-area of the square is 



^- = yi-^ = l + 0(A2). 

The area of the ellipse £ is irr'^, where r is the Q-distance between the fixed point (1/(2 
A), 1/(2 - A)) and the point of tangency Tq = ((1 + A)/2, 1). This is 



Thus, the area outside £ is 



A^\£ = l-!L + !Lx + oiX^). 
4 4 



Now we consider the dynamics, reverting to the representation (|17p. First, we calculate 
the area covered by the images of the regular atoms under the outer map L°'^^ . These 
atoms are polygons whose vertices are listed in appendix A.l. 



In the first-order area calculation, we can expand formulae (1521) and (I53p in Taylor 
series about A = 0, keeping only the lowest two terms and an O(A^) remainder. For 
m = 1, the quadrilateral A°"* reduces to a right triangle [(A + 2/3, 1), (1, 1), (1, 1 — A)] with 
area A/6-1- O(A^). Since the transit time under F is 7 (theorem [TT]) we get a contribution 
to the overall area sum of 

^A + 0(A2). 

For m = 2, the hexagon Ag*^* reduces to a quadrilateral with vertices {(| -|- A, 1), 
(| -|- A, 1), (1, 1 — A), (| -|- A, 1 — ■!)}, and hence area A/15. Multiplying by the transit time 
of 11, we get 

The remaining domains A^* are the quadrilaterals [PQ{m), Po{m—l),Pi(m—l), Pi{m)], 
L^'^^-symmetric about the respective axes [Po{m), Pi{m — 1)]. We use the formulae (jSOh 
and (ISTI) to calculate the area of A2"* as 



AT = \/r^AV4 (Po(m)-Pi(m)) X (Pi(m)-Pi(m-l)) 
A^cos((2m- 1)0) 
8nL-iSm((2m+(2A:-l))0)' 

Expanding the right-hand side in powers of 9, we have 

A^ {l + m'^e^r{e,m)) 



A out 

™ ~ 803(2^ _3)(2m-l)(2m + 1) 
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where, as we prove in appendix A. 8, r(6, m) is uniformly bounded above and below for 
3 < m < M, < < sin~^(A/2). Multiplying by the transit time 4m + 3 and summing 
over m gives 

^g(2m-3)(2m-l)(2m + l) 20 vr ^ ^' 

Em^{4m + 3) _ tt 

(2m-3)(2m- l)(2m + l) ~ 4A ^ ^ °^ ^' 



m=3 



and hence 



M(A) 

5](4m + 3Mr = ^A + 0(A2). 



20 

m=Ji 



Combining the results for all m, the area contribution is 



= Za + 0(A2). 
4 



Second, we compute the area covered by the orbits of the regular atoms in A^"^, under 
the inner map L"\ Here again we start with exact formulae for the vertices of the regular 
atoms given in appendix A.l, and the expansions of lemma [TH Each A™ is a quadrilateral 
[Qo{n),QQ{n — l),Qi{n — l),Qi{n)] with L°"* symmetry axis [Qo{n — l),Qi{n)] and the 
area is calculated to be [12] 



= - AV4 iQoin - 1) - Qo{n)) x (Qo(n) - Qi{n)) 

= ^t(1 + t2)A2 - 1(1 + 4t + 4r2 + At^ + 3r^) + rAA^ 
2 o 



where t] and r were defined in (164p . and is a function of A and r which satisfies < 29. 
Next we sum the atom areas multiplied by the transit times 4(n — 1) (theorem [6]) , and use 
the Euler-Maclaurin formula to turn it into an integral: 

N{X) 

= 5]4(n-lMjf 

n=l 

TT 

= lj'dv{^-X-r]) (4(t + t3) - (1 + 4r + 4t'' + 4t^ + 3t^)A) + ©(A^) 



Here we omit the detailed handling of error terms (Euler-Maclaurin endpoint and 
remainder terms, change of integration variable, adjustment of integration limits), which 
parallels that encountered in the calculation of disk areas in section [3j Details may be 
found in [12]. 

Now the regular atoms of L are the intersections of the 'in' and 'out' atoms — see 
equation (1591) . The mismatch between the union of the two sets of atoms is O(A^), from 
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theorems [6l and [TTl while the respective transit times are 0(A~^). Hence the area covered 
by the images of the regular atoms of L is ^™ + + O(A^). Combining the above area 
expressions gives us the desired identity 

□ 

Although the irregular atoms of L do contribute to the total area, the proof of theorem 
[23] avoids the great difficulty of deriving bounds for their return times. 



7 Extension of our results 

We consider, briefly and informally, some extensions of our results. 

The regularity of the limit A ^ 0^ depends on the very simple dynamics at A = 0, 
for which there is a single atom. There is one dominant fixed point, while all other cycles 
disappear at A = 0, and are arranged in a hierarchy that affords a perturbative analysis 
in the small parameter 9. Some asymptotic features depend on the product of 6 and the 
return time of relevant domains, so that large return times lead to an unavoidable non- 
uniformity in the convergence to the limit. However, this problem remains tractable, due 
to the regular arrangement of the atoms. 

Similar features are present in the limit A ^ 0~, so we expect similar phenomena. For 
— 1 < A < 0, there are three symbols: /-(x, y) € {0, 1, 2}. However, the atom VLq is a single 
point (the origin) and can be neglected. The atom Oi is again a quadrilateral with area 
0(1), while VI2 is a triangle with area 0(A), which plays the same role as the atom Q.q in 
the positive A case. 

For — 1 < A < 0, the basic cycles are 

t code denominator numerators radius 

1 (T) 2- A (1) xo(A) (87) 

2 (M) 4-A2 (4 + a, 2(1 + A)) l-xo(A) 

The fixed point is determined by the same rational function as for the A > case. The 2- 
cycle is different; however, changing the sign of A, together with a reflection with respect 
to the centre of J7, reproduces the same behaviour to first order in A. The symmetry 
properties of these cycles are also unchanged. 

The linked sector maps of the 1- and 2-cycles for A < are shown in figure [TUl 
Compared to the case A > 0, there is now a tighter connection between the two sectors, 
which results in shorter transit times. As for A > 0, one constructs the transit maps 
2^m,out^ and identifies a two-parameter family of regular fixed points with periods 

i(m, n) = 4(m + n) — 3 
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Figure 10: The linked sector maps of the 1-cycle and the 2-cycle for a small negative value of 
the parameter A. The corresponding sectors S and S' are represented in light and dark grey, 
respectively. Under the action of the map F'^, both sectors turns counterclockwise. The points 
in S that do not end up in S' (under F) comprise the triangle Ti; those that do comprise the 
quadrilateral T2. 

and symbolic codes 

^(n,m) ^ (l4n-l^(2,l)2'"-l) n > 1, m > 1, 

to be compared with ([60]) and ([6T]) . respectively. As before, there are anomalous cycles of 
total area 0(A), which interact with the above structure in regions of phase space of area 
O(A^). The main cycle of this kind is a fixed point located in the North- West corner of 
the square Q. 

Returning to the limit A ^ 0^, we recall that the mismatch between the sector maps of 
the 1 and 2-cycles generates a collection of small irregular domains, which were neglected 
in our perturbative analysis. The dynamics over these domains is determined by the 
interaction among certain periodic orbits lying outside the sector S. 

As an example, let us consider the turnstile associated to the domain Ti (figure [6]). 
Under the action of G o F~'^, the left triangle comprising the turnstile is mapped to the 
triangle $ = (0, 7, 8), of area O(A^), lying between the atom A°"* and the boundary of 
A (cf. the first equation in (f37j) ). 

The domain ^ belongs to the first regular atom of L™. The partition of $ into atoms 
under is shown in figure [11] for two parameter values approaching zero from above. 
The low-transit-time atoms (all symmetric) belong to three families, coloured blue, red, 
and green in the figure, and the picture suggests convergence to an asymptotic regime, 
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Figure 11: (a) The irregular domain $ on the eastern boundary of A^""*, showing three famiUes of 
atoms. Here X = (b) Same as (a), but with A = In both cases the horizontal and vertical 

dimensions have been rescaled to make the structure of the partition visible. 

although the convergence seems slow. The atoms' respective codes and transit times are: 



Colour 


Code 


Transit Time 


Blue 


(1(1101101010)^™-^ 110111) 


20m - 3 


Red 


(1(110)4^111) 


12m + 4 


Green 


(1(110)4™-^ (1101101010)(110)4™111) 


24m + 11 



We see that these codes shadow two periodic codes (110) and (1101101010). These cor- 
respond to a 3- and a 10-cycle lying outside S, whose cells have area O(A^) and 0{X'^), 
respectively. 

In addition to A = 0, there is trivial dynamics (a single atom) also for A = ±1 (rational 
rotation number p = 1/6, 1/3), and for these values we expect the near-rational dynamics 
to be broadly similar to that considered in this paper. By contrast, the dynamics for all 
other rational rotation numbers is highly non-trivial. The simplest cases correspond to 
quadratic irrational A- values (eight in all), for which the phase space is tiled by infinitely 
many periodic cells which admit an exact renormalization [23]. Preliminary investigations 
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for A = (\/5 — l)/2 (p = 1/5) indicate that the construct of linked sector maps can 
only serve as a local model for near-rational dynamics; the global properties result from 
complex interactions between different local models. We have observed a highly non- 
uniform convergence to the asympotic regime, due to the presence of unbounded return 
times in the limiting rational dynamics. This problem requires further investigation. 



Appendix 

A.l Rules of interval arithmetic with rational rounding 

Suppose we have rational lower and upper bounds on a set of real quantities ri, r2, . . . , r„ 
and we wish to obtain a similar bound on a rational function R of these quantities. We 
can do so by assigning to the quantity the closed interval [oj, hi] where ai < ri < bi are 
the given bounds. If r is a rational number, then we assign to it the closed interval [r, r]. 
We then proceed to deduce a bounding interval for R by performing, in any convenient 
order, a sequence of elementary arithmetic operations, taken from the following list: 

(z) [ai, 61] + [a2, 62] = [oi +a2,hi+ ^2], 

(m) [ai,6i] X [02,62] = [min{aia2, 0162, 6102, 6162}, max{aia2, 0162, 6102, &162}], 

{in) For a6 > 0, [a, 6]"^ = [min{i, i}, max{i, i}], 

[iv) For a > 0,6 > 0,n > or a < 0,6 > 0,n > 0, [0,6]'' = [min{a", 6"}, max{a", 6^"}], 

{v) For a > 0, 6 > 0, n > 0, [a, 6]-" = ([a, 6]")"^ 

The result is an exact, but obviously non-unique, pair of bounds for R. 

Since multiple application of these rules can lead to prohibitively large denominators for 
the resulting rational bounds, it is convenient to round off after every elementary operation. 
The lower and upper bounds are rounded down (resp. up) to the nearest rational number 
with a power of 10 in the denominator and a specified number of significant digits in the 
numerator. Any common factors in numerator and denominator can be divided out. For 
example, with the 4-digit rational rounding used in this article, the interval [7r/2, tt] would 
be rounded to [157/100,3142/10000]. 

The Mathematica interval functions used to obtain rigorous bounds throughout this 
article are listed in sec. E.l of [12]. 

A. 2 Proof of Pfc(m) estimates in lemma [141 

The initial terms in the stated formulae for PQ{m) and Pi{m) result from the Taylor 
expansion of (|50]) and (fSTI) up to second order in the x-component and third order in 
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the y-componcnt (in the latter case, the zcroth order term is trivially equal to unity). 
Our task is to establish the stated numerical bounds on the remainder coefficients ri, i = 
Ox, Ix, ly, Olx, lOy, lly. 

We begin with 

rocc = (2m + l)Po(m)^ - 1 - m - (2m + 1)A + ^(m + 1)X'^. (88) 
Defining auxiliary variables u and v by the expansion formulae 

2 4 
X X 

xcot{x) = 1- — - (l + u) — , 

Inserting these expansions reduces PQ{m)x, hence Vqx, to ratios of polynomials in X,m,u, 
and V. The numbers of terms is in these polynomials is enormous, but the expressions 
are easily manipulated with computer assistance. Moreover, all h\it a handful of terms 
have sufficiently high powers of A that they contribute insignificantly to the final estimates 
(recall our cut-off condition A < Aq = 10~^). The detailed calculations have been relegated 
to the electronic supplement [12] (E.7.1). 

To obtain a bounding interval for tqx, we use the following bounds [12] (E.7.1) for our 
variables: 

A G [0, 10"''], m G [2, \~'-^% u G [-3/2500, 3/1000], G [0, 2 x 10"% 

Each monomial of the form ■iz\""m}'u'^v^ is assigned a single bounding interval using the 
rules of interval arithmetic (see appendix A.l) with fixed precision rational rounding down 
(resp. up) of lower (resp. upper) bounds. The monomial bounds are then be rigorously 
combined, again using interval arithmetic, to give the final bounds for tqx stated in the 
lemma. The same technique was used to establish [12] (E.7.1-2) the stated bounds on the 
other remainder factors in the lemma. 

We have implemented the process described above with the aid of Mathematica func- 
tions involving only elementary logical and arithmetic manipulations of symbols and inte- 
gers. The specific procedures are listed in the Electronic Supplement [12], sec.E.l. Note 
that while our estimates are rigorous, they are not optimal, and can easily be changed by 
permuting the order in which monomial bounds are combined. 

A potential difficulty of the method comes from the fact that the upper bound for m 
grows as A~^/^, which is numerically unbounded for A in the chosen range. The problem 
is avoided by a simple trick based on the fact that all polynomial expressions which we 
wish to bound have only monomials of the form ztA'^m'' • • • where b is either negative or 

no greater than 4a. In the first case we use m'' G [0, 2^*] , and in the second case we use 
A'^m^ G [0,Ao"''/^]. If the original rational function is not in this form, we try, if possible, 
to make it acceptable by dividing both numerator and denominator by a common power 
of m. 
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A. 3 Proof of Qfc(n) estimates in lemma [T4l 



Taylor expansion of ()54p and (j55p has been carried out with respect to A with r held fixed. 
The remainders, expressed as rational functions of A and r are 

r2x = = rsy = (89) 



For r2x, we have 

AA2^ = -64 + 64r2 + 32A + 32tA - 32r2A - 32r3A + 32A2 - IGrA^ - 40t'^\'^+ 

IGr^A^ + Sr^A^ - 24rA3 + IGr^A^ + Sr^A^ - 4rA^ + lOr^A^ - 4t3A^ - 2r4A^+ 
V4- A2(32 - 32r2 - 16A - 16rA + IGr^A + IGr^A - 12A2 + 8t\'^+ 
16r2A2 - 8r3A2 - 4r4A2 - 2A3 + SrA^ - Gr^A^ + 2t^\^ + r^A^), 

P2x= 32A2 (^-8r + 2A2r + ^/4^^(A- Ar^)) . 



Application of interval bounds to these expressions runs into a serious difficulty: the 
denominator has a zero (barely!) within the allowed range < r < 1, namely at 



2V4^-4-A2 a ^,^3, 

Fortunately, the numerator shares the same zero, and we can divide both J\f2x and T>2x by 
the factor (r — ti). The result takes the form 

r2x = ^ (2V4-A2 - 4 + A^) , (90) 

where I?2x ^^'^ (especially) are lengthy polynomials in A, r, and \/4 — A^. 

To apply the interval methods used successfully in appendix A. 2, it turns out that we 
need to use a high-order estimate for the square root, namely 

A2 A^ A6 5A8 7A10 21vX'^ 



^4 - A2 = 2 



4 64 512 16384 131072 2097152 
Inserting the interval estimates 

1117 

""^[^'1000^' ^^[O'l]' ^e[o,Ao], 

in each of the 3 polynomial components of r2x in (|90p and applying the rules of appendix 
A.l to combine them, we arrive, finally at the stated bounds. Applying the same tech- 
niques to r^x, r^y, r4x, r^x, we get the remaining estimates in lemma [TH Details of the 
calculations may be found in [12], sec. E.7.3. 
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A. 4 r estimates 



In addition to the contents of lemma O concerning the cartesian coordinates of the A'^^ 
vertices, we will also need information relating to their r-coordinates (see (|64p ). Estimates 
for the latter may be obtained by combining the formulae of lemma [HI In the follow- 
ing lemma, the remainder terms are rigorously bounded [12] using the same techniques 
discussed in appendices A. 2 and A. 3. 

Lemma 24 For A < Aq = 10~^ and 2 < m < X~^/^, we have the expansions 

TiPoim)) = + ^ (5 + 2^-1 + ^-2) 

+— (24 + (57 - 16m(m + 1))^"^ + 12^-^ + 3fi-^) + ro^A^ (91) 
4o 

where /i = 2m + 1, ror G [-fisS' TSm] > "■^'^ 

+— (24 + (129 - 16(m + 2)(m + 3))//'"^ + eO/u'-^ + 75//'-^) + rir\^2) 
48 

where /u' = 2m — 1, rir G [~^) ^] • Moreover, 



r 383 2819 ] 
L25OO' 1000 J 



with tqit- G 

As a corollary, we get bounds on T{Qo{n)) restricted to a single value of m. 
Corollary 25 For A G [0,10-'^],m G [3,X-^/'^], 



mT{Qo{n)) = rr £ 



2 1001 



5 1000 



Proof. For those atoms A™ which intersect A^*, we have the inequalities 

r(Pi(m - 1)) > r(Qo(n)) = tan(7(go(n - 1)) - 26) > (93) 
where, for A < 10~^, 

tan(20) = ^^5^ < (1 + 10^«)A. 
2 — A 

From lemma [Ml it is not difficult to derive [12] 

. 1 /5 \ , 1467 

T{Po{m)) = — + - + ro A, ro > 



2m + l ' l4 ' V 2 X 106' 
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T(Pi(m)) = + - + r2 A, rs < 



2m -1 \A y ' 500 

With the inequahties (vahd in the cut-off regime) 

m 3 m ,1 
> < 1, 0<Am< 



2m + l~7' 2m-3~ ' ~ "1000 

we obtain 

- < mT[Qo{n)) < 



5 ' ^ " 1000 

□ 



A. 5 Proof of lemma [161 



The differences in the n-coordinate are related to those in r by 

n2 -ni = -i(tan-i(ri) -tan-i(T2)) = . _^ tan"^ (t^^) ■ (9^) 
2t) 2sin (A/2) \^ + 'TiT2/ 

From lemma [2^ we have 

r(Pi(m)) - r(PoM) = ^^^^ + 0(A), 

r(A(m - D) - riMm - 1)) = _L + 3 + ^^^^ 

^A<r(Po(m-l))-r(Pi(m))<l^A. 

Inserting these estimates in (j94p and considering the small- argument behaviour of the 
inverse trigonometric functions, we obtain the desired relations. 



A. 6 Proof of lemma [TH] 



The method of proof is a familiar one. The remainder factors rzx and rzy are constructed 
by combining (|69p and (j70|) with the insertion of PQ{m), Pi{m — l),QQ(n — 1), and Qi{n) 
from lemma [TSj Here care must be taken to express Qo{n — 1) in terms of the variable 
r = T{Qo{n)): 

r{Qo{n - 1)) = tan - InO + 29^-^ + ^""^'^^ 



2 7 1-Ttan(26')' 

The resulting expressions (see [12] for details) for rzx and rzy are enormous polynomials in 
A, m, r, and the various remainder factors appearing in the expansion formulae for the atom 
vertices. Because the Z(m, n) lies within a single atom A^^, we improve our estimates by 
replacing r by r/m, where, from appendix A. 5, r £ [2/5, 1001/1000]. Application of the 
rules of interval arithmetic (appendix A.l) then yields the uniform bounds stated in the 
lemma. 
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A. 7 Crossover calculation 



We now return to the crossover phenomenon for a more precise treatment. First of all, 
we observe that the angle between the line segments [PQ{m), Pi{m)] and [QQ{n{PQ{m)), 
Qi{n{PQ{m))] changes sign as one proceeds from right to left in A. For some m = m* ,n = 
n* = n{PQ{m*)), the segments coincide and the atoms A'^i and A™ have their vertical 
edges nearly parallel. We adopt m = m* as our precise definition of the crossover. By the 
construction of section [3l the orientation of the line through [PQ{m), Pi{m)] differs from 
the vertical by a generalized rotation C"^™", while [Qo{iT')iQi(''^)] is rotated by C^"^^. 
Thus crossover corresponds to 

^4r7i*+4n*+3 

Obviously there are infinitely many solutions, but the relevant one, corresponding to a 
single crossover as m varies from 1 to Af, has a negative sign, with m* satisfying 

(4m* + 4n{Po{m*)) + 3) sin" 

where n{Po{m*)) is calculated by solving Qo{n) = X for n, then setting X = PQ{m*). 
The equation for m* can be solved perturbatively to obtain 

^*(A) = ^fi + ^ + il + ^ + 0(A^)V (95) 
V2\ V 3 120 1008 ^ ^ ' ^ ^ 

■*(\\ ^ 1 3 ^/A vrA A3/2 

= 4A-71X"4-Vl"^"TWl + ''^'^ 
In particular, m*(A) diverges to infinity as A — > 0. 




A. 8 Proof of area estimate 



m 



We are seeking uniform bounds, for 3 < m < M and < ^ < 10 ^ on 

e-'^m''^c{{2m - 1)9) 
r{e,m) - ^^^2^ _ 3)^^g^^2m - l)e)s{{2m + 1)6) ~ ' 

where 

c(x) = cos(x), s(x) = . 

X 

Note that in calculation of bounds we will not be able to benefit from the previous cut-off 
on m. 

Following our usual practice, we define auxiliary variables u and v via 

, , , (1 — ti)x^ , , , (1 — v)x^ 

c(x) = l ^- — , s(x) = l ^- . 

^ ^ 2 24 ^ ^ 6 120 

Bounds < u,v < 1002/1000 are straightforward to establish [12] (sec. E.11.1). Insertion 
of these expressions reduces r{6, m) to a ratio of polynomials in the variables 9, m, u, and 
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V. Although each polynomial can be uniformly bounded above and below using the same 
methods used throughout this article, there is a serious problem: the denominator range 
includes zero. On the other hand, numerical evaluation of the denominator shows that it 
is almost certainly negative definite. 

To proceed, we isolate the leading contribution consisting of terms proportional to 
m^O^ (for all others, the power of m is less than that of 0). There is no u dependence of 
these terms. With respect to v, we make the worst-case assumption, setting it equal to 
zero for all negative terms and to its maximum value of 1002/1000. The leading part of 
the denominator then reduces to 

Vi{e,m) = h{w?e'^), 



hix) = -1728000 + 3456000X - ^^^i^^^x^ + 1433600a; 

5 

2328576 4 384768768 5 24096096064 g 
25 ^ 3125 ^ 1953125 ^ ' 



,r2 



We need to show that this polynomial is negative definite over the range < x < ^ < 
xi = iQQQQ ■ Direct application of interval arithmetic fails, as expected, producing a range 
of values containing zero. On the other hand, the same method, applied to the derivative 
h'{x) produces a pair of positive bounds.Thus h{x) is a monotone increasing function. 
Since it is negative at both and x\ it must be negative everywhere on the interval. 
Specifically, 

Vi{e,m) G [-1728000,-110000] 

Meanwhile, the numerator N{9, m, u, v) and the non-leading part of the denominator, 
V2{6,m,u,v) can be bounded by our original method, obtaining [12] (sec. E.11.1) 

J\fi9,m,u,v) £ [-24730000,27090000], V2i{e, m,u,v) G [-31800000, 14920]. 



Combining the bounding intervals according to the rules of appendix A.l, we get, 
finally 

N r 1301' 
rie,m)£ -285,—- 
5 



□ 
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